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STRATIGRAPHY AND RADIOCARBON DATES 
AT SEARLES LAKE, CALIFORNIA 
RICHARD FOSTER FLINT and W. A. GALE 


LBSTRACT. Upper units of the strata beneath Searles dry lake, in a desert basin in 
outheastern California, record two deep lakes each followed by precipitation of evaporites 
d at least near-desiccation, Each lake is thought to imply a pluvial climate; the times 
near-desiccation represent the warmer, drier, nonpluvial climate of today. C* dates of 
amples from these strata show that the later of the two pluvial climates lasted from before 
23,000 to ~ 10,000 yr. B.P.; it was therefore contemporary with the classical Wisconsin 
glaciation of central United States. The earlier pluvial had been well established by 
~ 46,000 B.P. and had begun to wane before ~32,000 B.P. The later history of the lake 
n the Searles basin is in agreement with the published data on the later history of Lakes 
Bonneville and Lahontan. Still earlier lake fluctuation can be discerned in underlying 
trata, but these are not controlled by dates. 


LOCATION AND BRIEF DESCRIPTION 

Searles Lake occupies one of the many structural basins in the Mojave 
Desert region of southeastern California; it lies 170 miles northeast of Los 
Angeles (fig. 1). It is notable for the great volume and purity of the salt beds 
occurring beneath its floor. : 

Like many other “lakes” in the Basin-and-Range province, Searles is essen- 
tially a dry lake (plate 1). It has an area of about 60 square miles, four-fifths 
of which is salt-incrusted sand, silt and clay. The central part, about 12 square 
miles in area, has a firm salt surface from which common salt can be harvested. 
The salt body, the surface of which is exposed in its central part, is thick and 
grades radially outward, beneath peripheral accumulations of mud, into salifer- 
ous clay, silt and sand, and altogether underlies an area of about 40 square 
miles (fig. 2). This enormous mass of salines is porous; its voids, estimated at 
0) percent of its volume, are filled to the surface with a saturated bittern, Al- 
though the bittern varies in composition, its chief constituents are sodium and 
potassium chlorides, sulfates, carbonates, and borates. Following infrequent 
floods, pools on the surface of the salt are colored red owing to the abundant 
presence of red bacteria and algae. 

The bittern is worked on a large scale by two companies, the American 
Potash & Chemical Corporation at Trona on the northwestern margin of the 
lake, and the West End Chemical Company at Westend on the southwestern 
side. The former plant produces potash salts, borax, boric acid, soda ash, salt 
cake, bromine, lithium carbonate, and phosphoric acid, while the latter firm 
produces mainly soda ash, borax, and salt cake from the bittern. 

The purpose of the present paper is to describe in general terms the 
stratigraphy of the sediments beneath Searles Lake, to present C** dates ob- 
tained from core samples, and to comment on the significance of both strati- 
graphy and dates for the Pleistocene history of the Searles basin. 
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Fig. 1. Map and cross-section showing Searles Lake and chain of pluvial lakes 
(shaded) of which Ancient Searles lake was a part. Dotted lines show routes of former 
connecting streams. Sources: H. S. Gale 1915; Blackwelder 1933; U. S. Geol. Survey 
topographic maps. 


THE CHAIN OF BASINS AND ANCIENT LAKES 

The former existence of a large lake in the Searles basin is established by 
strandlines, subsurface stratigraphy, fossil fish and mollusks, and deposits of 
tufa. The former lake was recognized by H. S. Gale (1914a; 1915) as one of a 
chain of four lakes—the ancient Owens, Searles, Panamint, and Manly lakes 
(fig. 1), on which other published references are Blackwelder (1933; 1941), 
and Antevs (1938). The chain was fed principally at its northwestern end, by 
streams draining a 130-mile segment of the east slope of the Sierra Nevada. 
The ancient lakes have been generally regarded as dating from episodes of 
glaciation, when lowered temperatures reduced evaporation, and increased 
discharges resulted from increased precipitation and secondarily from the 
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presence of numerous glaciers on the east side of the High Sierra. As a result 
the lakes filled their basins and overflowed from one basin to the next. 

The outflow from Ancient Searles lake discharged into the Panamint 
Valley, where another lake, having an area of 272 square miles and a maxi- 
mum depth of more than 900 feet (Thompson, 1929, p. 186), discharged east- 
ward through Wingate Pass into Lake Manly in Death Valley, the terminus of 
the lake system. Lake Manly, described by Blackwelder (1933), had a maxi- 
mum depth of about 600 feet, but was relatively short lived. As it had no out- 
let, it fluctuated at all times with varying inflow and evaporation, Probably it 
was nearly fresh during its expansion and became increasingly saline during 
its shrinkage. Because of the state of preservation of strandline forms and de- 
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Fig. 2. Map of Searles Lake showing section lines, approximate areas of exposed 


and buried salt, and locations of 30 toreholes (from among several hundred), mentioned 
in the text or shown in figures 3 and 4. 
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posits, and of the streamway at Wingate Pass, Blackwelder inferred that Lake 
Manly and Ancient Panamint lake were contemporaneous with the highest 
shore features in the Searles basin. but not with lower ones. In other words. 
that with cessation of outflow of Ancient Searles lake the chain consisted only 
of the Owens-Searles system, although of course small independent fluctuating 
lakes might have existed then in Panamint and Death Valleys. The chain of 
lakes, as described and interpreted in the literature, affords a reasonable ex- 
planation of the stratigraphic section, discussed hereafter, beneath the floor of 
Searles Lake. 

All four basins contain successions of strandlines, but those in the Searles 
basin’ are probably more numerous and more distinct than those in the other 
basins. They include narrow wave-cut benches and pebble beaches, so dis- 
continuous that no single one continues around the basin without a break. The 
highest and strongest of the conspicuous strandlines is at an altitude of about 
2,260 feet. On the southeast side of the basin it leads into a pass that evidently 
was the lake outlet while the water surface stood at that altitude. 

Blackwelder (1941; 1954, p. 35-36) reported from reconnaissance the 
possible presence of two groups of strandlines: a weathered and dissected group 
that he thought might correlate with the Tahoe glaciation in the Sierra Nevada. 
and a lower, fresher group that might correlate with the Tioga glaciation, 

When at its maximum recorded height the Ancient Searles lake had an 
area of about 385 square miles and probably was at least 750 feet deep. It was 
common to three basins, from west to east: Indian Wells Valley, Salt Wells 
Valley, and Searles basin proper (fig. 1). Indian Wells Valley is much shallow- 
er than Searles basin and contains a large volume of sand and silt, interpreted 
by H. S. Gale (1915, p. 269) as a deltaic deposit built by the inflowing stream. 
The present writers agree with Gale. and believe that he did not sufficiently, 
emphasize the occurrence of lacustrine silt in the southwestern part of Searles 
basin. The large body of sandy sediment in Salt Wells Valley continues east- 
ward with decreasing grainsize and altitude, partly fills the broad depression 
connecting Salt Wells Valley with Searles basin, and underlies dissected benches 
throughout an area aggregating many square miles in southwestern Searles 
basin. In that area it consists predominantly of dusky yellow green silt and 
fine sand resembling the coarsest fraction of the Parting Mud, described here- 
after, beneath the central part of the basin. Although it seems likely that the 
silt passes eastward beneath the floor of the basin and is continuous with the 
Parting Mud, we do not know that this is the case. 

The silt contains mollusk shells, Small collections of the shells from surface 
exposures in Owens basin, Indian Wells and Salt Wells Valleys, and Searles 
basin were examined by Dwight W. Taylor (written communication). He 
found that the nearest occurrence alive of some of the 13 species represented 
is in the Sacramento Valley in northern California, and on this basis inferred 
that when the fossil fauna was living the Searles Lake region had considerably 
cooler summers than it has today. The difference in mean July temperature 
between the two districts today is about 4.5°C. 


* Possibly first described by E. E. Free, in unpublished notes quoted verbatim by Hunting 
ton (1914, p. 559-560). 
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Fig. 3. North-south section AA’ (fig. 2), constructed from core-drill logs to show 
stratigraphic units described in text. 
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STRATIGRAPHY 

Searles Lake is underlain by a remarkable sequence of sediments, from 
which much has been learned about the Pleistocene history of the basin and 
from which further information will surely be derived. Even from the few 
borings available to him, H. S. Gale (1914a, p. 404) recognized that the sedi- 
ments are of lake origin. Since the time of his study many core borings have 
been made by the American Potash & Chemical Corporation and its predeces- 
sors, by the U. S. Geological Survey, and by others. Numerous borings have 
been made to depths of 70-90 feet, many to about 150 feet, several have pene- 
trated to approximately 300 feet, and one boring extends to 875 feet. The lo- 
cations of 30 boreholes are shown for example in figure 2; two sections con- 
structed from the logs of eighteen of these borings are shown as figures 3 and 4. 

Study of the resulting cores has confirmed H. S$, Gale’s opinion and has 
made it clear that the Searles Lake stratigraphy represents lakes that succeeded 
each other in time and that fluctuated in response to changes of climate. It is 
our purpose not to detail the stratigraphy but to describe major units. Those 
which are of chief present interest, because they are young enough to be 
datable by C'*, occur within the uppermost 200 feet of the section, and are 
shown by a graphic log, figure 5. The names used in this paper are those used 
informally by the American Potash & Chemical Corporation. These are, from 
the top downward, Upper Salt, Parting Mud, Lower Salt, Bottom Mud, and 
“Mixed” layer. The descriptions and thicknesses given below pertain to the 
central part of the basin. The units change facies toward the margins, but we 
are not concerned here with their marginal relations. As the stratigraphic units 
are believed to include two natural pairs, a Parting Mud-Upper Salt pair and 
a Bottom Mud-Lower Salt pair, they are so described in the following para- 
graphs. 
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Fig. 4. West-east section BB’ (fig. 2). 
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PARTING MUD AND UPPER SALT 
The Parting Mud, 10 to 13 feet thick, consists of saliferous clay and fine 
silt, dusky yellow green when dry, with parallel, flat-lying laminae (pl. 2, fig. 
PLATE 1 
Searles Basin and Lake 


American Potash & Chemical Corporation 

Fig. 1. Searles basin seen looking southeast from directly over the Argus Range, the 

western rim of the basin, The distance between Trona and Westend is about 4 miles, Out- 
let of the ancient lake is out of view to the right. 


Ward C. Smith, U. S. Geological Survey 


Fig. 2. Shoreline of Ancient Searles lake along east side of basin. 
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1). Seen megascopically the laminae consist of alternating light-colored layers 
mostly less than 2 mm thick, and thicker dark layers up to 1 cm thick, Under 
90x magnification the laminae show no evidence of possessing graded bedding; 
instead the darker laminae are seen to consist of distinct laminations of a small- 
er order of size, some of which are lenticular and can not be followed laterally 
through the length of a hand specimen. Examination, by Dr. J. E. Sanders, 
Yale University, of a single thin section under a petrographic microscope re- 
vealed that the small coarse fraction, fine sand to coarse silt in size, consists of 
quartz with small amounts of biotite, tourmaline, glass, and unidentified 
minerals. The finer fraction, a large proportion of the whole, includes abundant 
unidentified carbonates. The proportion of clay-size to coarser particles was 
not determined. The dark-colored layers include significant amounts of clastic 
material; the light-colored layers appear to be chemical precipitates. The 
samples from which two C' dates were obtained (table 2, GRO-1460, 1464) 
contained only about 12 percent insoluble, presumably clastic material (HI. de 
Vries, unpublished), most or all of which probably resided in the dark-colored 
layers. Whether the substances in the light-colored layers were sedimented in 
part from mechanical suspension or were all precipitated chemically on the 
floor of the lake was not determined. In either event they appear to have un- 
dergone some recrystallization. 

Although the material, here called mud, itself consists predominantly of 
salts, its content of true evaporites is not known. To the extent that the sub- 
stances in the light-colored layers are chemical precipitates on the lake floor, it 
can be inferred that the lake water was saturated with carbonates even during 
times when the mud was accumulating. There is little doubt that two substances, 
at least, in the Parting Mud are true evaporites: (1) disseminated crystals of 
the relatively insoluble minerals gaylussite and pirssonite, particularly near the 
upper and lower contacts, and (2) large crystals of borax (tincal) scattered 
sparsely through the uppermost few inches of the mud, immediately below the 
layer of trona that constitutes the basal part of the Upper Salt. Probably the 
borax was formed through seasonal chilling of the concentrating brine just 
before saturation with trona was reached. 

Apart from the occurrence of these crystals, the top and base (plate 2, fig. 
2) of the Parting Mud are in sharp contact with the overlying and underlying 
salt bodies. 

The Parting Mud contains fossil fishes, pollen, and disseminated organic 
matter. It was examined* for diatoms, but with negative results, presumably 
because during diagenesis the strongly alkaline water dissolved any diatoms 
then present. The Mud will also be examined,’ for limnic animals in the hope 
that data on salinity, temperature and depth can be derived from them. 

A fossil fish 8 inches long was found at a depth of 68 feet in the Parting 
Mud in core ff (fig. 3). It was identified* as a member of the minnow family 
Cyprinidae, probably of the genus Siphateles. No significant ecologic inference 
has been drawn from this find. 

Pollen from both the Parting Mud and the Upper Salt has been examined 
? By Dr. Ruth Patrick, Academy of Natural Sciences, Philadelphia. 

* By Dr. E. S. Deevey, Yale University. 


* By Dr. Carl L. Hubbs, Scripps Institution of Oceanography. 
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in detail (Roosma, 1958). The results show, through the stratigraphic distribu- 
tion and relative abundance of juniper compared with artemisia and chenopods, 
that the local climate was cooler and wetter during deposition of the Parting 
Mud than during precipitation of the overlying Upper Salt. 

Chemists of the Research Laboratory of the American Potash & Chemical 
Corporation, under the direction of W. A. Gale, found that the Parting Mud 
contains disseminated organic matter varying in amount from less than 0.1 
percent to as much as 1.5 percent of the weight of the sediment. Samples were 
concentrated for study in this way: sections of drill core from selected depths 
were broken up, dried at 140°F, and crushed to a coarse powder. Organic 
material was extracted in a large Soxhlet-type extractor by dissolving in 
methanol. The clear, filtered methanol solution was evaporated down to a 
syrup, and distilled water added drop-wise over a period of several days to 
precipitate the organic matter in particles coarse enough to be separated by 
filtration. The residue was washed with dilute HCl and then with distilled 
water. The precipitate was dried and dissolved in acetone, and the small in- 
organic residue was filtered off. The acetone filtrate was evaporated down to a 
thick slurry on a water bath and further evaporated to dryness in a vacuum 
desiccator. From the final product. weighing 50 to 70 gm, a small sample was 
analyzed for ash content by decomposing with sulfuric acid, igniting at 800°C, 
and weighing the residue. In each sample the residue contained molybdenum, 
iron, and phosphate, but no alkali metal salts such as would have been present 
if soluble inorganic carbonates remained in the organic extract. 

The organic matter in dry form is almost black in color, but in dilute 
solvent solution it is an orange dye. It possesses pronounced fungicidal prop- 
erties. Because of its color and its molybdenum content it is believed to have 
originated, in part at least, in some of the various forms of chromogenic or- 
ganisms such as the red and green bacteria now living on the surface of 
Searles Lake. These red bacteria have been found to contain comparable 
amounts of molybdenum. Such acetone-soluble organic matter constitutes, on 
an average, about 0.2 percent (by weight) of the Parting Mud, with some thin 
layers containing more than 1 percent dry basis. 

The organic matter is a complex mixture of substances. An investigation 
(Vallentyne, 1957) using chromatographic methods has shown the presence 
of a number of carotene compounds, xanthophylls and possibly echinenone, 
together with several other unidentified substances. The Parting Mud was also 
found to contain green chlorophyll-like pigments, extractable with methanol 
but not with acetone. Hence the organic content of the sediment seems to be a 
composite residue of algae, bacteria, and other micro-organisms that infested 
the waters of Searles Lake during deposition of the Parting Mud. 

The Upper Salt, 60 to 90 feet thick, consists of layers of salts, mostly 
alkali chlorides, carbonates and sulfates. Most of the salt layers are rather 
pure, but some are separated from others by partings of mud. Evidently the 


Upper Salt consists of the evaporites precipitated during progressive desicca- 
tion of the deep lake in which the Parting Mud had been deposited. 

The basal part of the Upper Salt, immediately above the Parting Mud, 
consists of a layer of trona about 5 feet thick. As trona, the stable form of 
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sodium carbonate in equilibrium with atmospheric CO., has a rather low 
solubility in sodium chloride solutions, it would have been one of the first 
sodium salts (other than the calcium-sodium carbonate double salts such as 
gaylussite) to precipitate from these brines. The other evaporites are difficult 
to explain in physical-chemical terms, owing to the complexity of the sequence. 
Chatard (1890, p. 59; see also H. S. Gale, 1915, p. 286) experimented 
with the evaporation of water from Owens Lake under natural open-air condi- 
tions. Chatard’s analysis of the resulting evaporites, at various stages of the 
process, showed that (1) sodium bicarbonate was precipitated early; (2) car- 
bonates as a group were precipitated almost continuously throughout much of 
the process; (3) sulfate precipitation tended to concentrate during an inter- 
mediate phase; (4) sodium chloride precipitation increased progressively with 
time; (5) carbonate could be precipitated by chilling. The sequence of minerals 
in the Upper Salt is generally compatible with Chatard’s results, if it is as- 
sumed that initial compositions were comparable. Variations in the sequence 
admit the possibility that variations in climate occurred during the process, but 
are not known to demand them. 

As a pair the Parting Mud and Upper Salt are the record of two phases 
of evaporation of a deep lake. The rather sharp contact of the Parting Mud 
on the underlying Lower Salt (mostly halite and trona near its top) suggests 
that the deep lake formed rather suddenly, and rapidly became deep enough to 
permit the sedimentation of clay on the floor, in parallel, nearly undisturbed 
laminae. The cores examined afford no good evidence that the surface of the 
underlying salt was dissolved appreciably by the lake water overlying it. If no 
solution did in fact occur, this might be explained by the supposition that the 
basal laminae of the Parting Mud were deposited rapidly enough to form a 
quick seal over the salt. Even in this early phase evaporation must have oc- 
curred to an extent sufficient at times to saturate the water with relatively in- 
soluble pirssonite and gaylussite. The greater rarity of these minerals in the 
stratigraphic middle portion of the Parting Mud suggests that net evapora- 
tion reached a minimum around a midpoint in the history of the “freshwater” 
lake. 

Little information on salinity is obtained from the presence of Siphateles 
)} in the upper middle portion of the Parting Mud, because the members of ~~ 


(? 
this genus live under a wide variety of physical and chemical conditions in the 
region today (Miller, 1949, p. 451). 

It is not known whether the lake, during its Parting Mud phase, over- 
flowed its basin or not, because there is no physical connection between bottom 
sediments and strandlines. The stratigraphy can be reasonably explained on 
either assumption. 


That net evaporation increased towards the end of Parting Mud time is 
suggested by the increase in scattered pirssonite and gaylussite crystals in the 
top part of the mud. The sharp contact of Upper Salt on Parting Mud records 
the cessation of deposition of mechanically transported sediment, followed by 
increase in salinity of the water to the saturation point for trona. As the 
principal if not the sole source of sediment in the center of the Searles basin 
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was the Owens River-Owens Lake system, it is reasonable to suppose that the 
Parting Mud had been derived from that system, and that the clay fraction at 
least had remained in suspension while passing through Owens Lake. A drill 
core from the dry floor of that lake (Smith and Pratt, 1957, p. 5) showing 
more than 900 feet of clay and other mechanical sediments with virtually no 
evaporites, is compatible with this supposition, which. however. can be verified 
only through detailed study of the sediments beneath both basins. Very likely. 
then. the top of the Parting Mud records the cessation of overflow of Owens 
Lake. cutting off the chief source of both water and suspended sediment for the 
Searles basin. Cessation could have resulted most reasonably from warming 
and drying of the regional climate. a trend indicated by pollen study. whereby 
the change from juniper, prominent throughout Parting-Mud time. to desert 
plants, prominent throughout Upper Salt time, is shown to have been rathet 
rapid (Roosma, 1958). The indicated climatic difference is compatible also 
with the cooler, wetter climate implied by the fauna of the silt exposed in the 
southwestern part of the Searles basin, on the assumption that that silt is con- 
temporaneous with the Parting Mud. 

The cessation of clay sedimentation and the inception of predominant 
precipitation of trona were two separate events, although both are logically 
attributable to a single cause, cessation of inflow. With consequent marked in- 
crease in net evaporation, saturation with trona occurred in Ancient Searles 
lake, but the physical evidence affords no information on the time lag between 
the two events. Nor are the depths of water at the times of these events known. 
because the strandlines have not been correlated with the floor stratigraphy. 

The subsequent history of the lake is one of continued evaporation. and 
precipitation of evaporites. with changes and interruptions as suggested. al- 
though not clearly defined, by the sequence of substances in the Upper Salt. 
The thin layers of mud occurring within the Upper Salt are continuous enough 
through the central part of the basin to suggest the possibility of temporarily 
renewed inflows of sediment-laden water. Whether the mud in these layers 
originated upstream from Owens Lake and arrived via that lake. or whether. 
like the surface mud that overlies the periphery of the Upper Salt. it was 
washed in from the sides of the Searles basin itself. could be learned only from 
detailed examination of the mud. The possibility of renewed overflow of Owens 
Lake is discussed in a later section. in the light of C'* dates. 

In connection with inflow to Searles basin, calculation suggests that an 
influent stream with only modest discharge would suffice to balance evapora- 
tion from the surface of a lake in the Searles basin. According to Blaney 
(1957) about 79 inches of water can be expected to evaporate annually. undet 
the existing climate, from the basin of Silver dry lake, 75 miles southeast of 
the Searles basin. The climates over the two basins do not differ greatly, A lake 
occupying the entire floor of the Searles basin would have an area of about 60 
square miles, Assuming this area and the stated rate of evaporation. such a lake 
would lose about 11 billion cubic feet of water annually, To balance this loss a 


continuous influent discharge of only about 325 cu.ft./sec. would be required. 


Under a cooler, moister climate evaporation would be less and the inflow 


needed to balance eit would be red Tt ed correspondingly. 
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BOTTOM MUD AND LOWER SALT 

The Bottom Mud and the Lower Salt constitute a related pair similar in 
many respects to the Parting Mud-Upper Salt pair. Together they represent an 
older, relatively deep lake that subsequently evaporated to a shallow, very 
saline lake or perhaps to complete dryness. 

The Lower Salt, 25 to 35 feet thick, consists of layers of mixed salts with 
intervening partings of mud. Except that it is thinner and includes more mud 
partings, it generally resembles the Upper Salt, with trona and other carbonates 
predominating in its lower part and with halite occurring in its upper part, Its 
lower and upper surfaces of contact are sharply defined. The sharp lower sur- 
face is believed to represent the cessation of deposition of fine suspended sedi- 
ment and hence probably a time when Owens Lake ceased to overflow, The 
distinct upper surface resembles that of Searles dry lake today, and may there- 
fore represent complete desiccation; on the other hand it may represent only 
the floor of a very shallow saline lake. 

The Bottom Mud, like the Parting Mud, is overlain abruptly by the trona 
in the base of the succeeding salt body. Its base, however, differs from the 
clearly defined base of the Parting Mud in that the mud grades downward, 
through alternations of mud and evaporites, into the underlying “Mixed” 
layer. For the purpose of this discussion the base of the Bottom Mud is fixed 
arbitrarily at depth 219 feet in core L-W-D, following the usage of G. I. Smith 
(unpublished), the depth at which evaporites become more abundant than 
mud in the downward direction. As thus defined the Bottom Mud is about 100 
feet thick in the central part of the basin. 

From inspec tion of core Jogs this layer is judged to be essentially similar 
to the Parting Mud, although it is somewhat more variable in composition. A 
noteworthy difference is a layer of mirabilite (Na,SO,-10H.O), 5 to 11 inches 
thick, at about the stratigraphic middle of the Bottom Mud. The mirabilite 
implies marked chilling of the lake water, and is compatible with the concept 
that the deep lake resulted from climatic cooling that reached its greatest in- 
tensity near the middle of Bottom-Mud time. Whether such cooling did in fact 
occur depends partly on evidence that may be found in other areas. 


“MIXED LAYER 
Phe “Mixed” layer, gradationally underlying the Lower Salt, extends 
down below the bottom of the deepest hole drilled (borehole L-W-D: Smith 
and Pratt, 1957, p. 25-51); it is therefore at least 700 feet thick. In the upper- 
most 150 feet (depth 177.4 to 310-315 in core L-W-D) it consists mainly of 
relatively insoluble minerals such as nahcolite, trona, and gaylussite, with a 


little clay and no halite. Apparently this sequence represents a time during 
which the lake was saline but did not reach the crystallization point for NaCl. 
This condition could be explained by the assumption that at the time the lake 
was overflowing continuously into Panamint Valley. Conceivably the halite en- 
countered beneath the floor of Panamint Valley (Smith and Pratt, 1957, p. 58) 
at depths exceeding 80 feet represents solutes that entered from Searles Lake 
at this time. 


bd 
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The uppermost 150 feet of the “Mixed” layer are underlain by a sequence, 
more than 500 feet thick, in which trona and other carbonates are predominant, 
with much smaller amounts of halite and some layers of mud. Most of the 
zones are thin, but one 20-foot zone and one about 70 feet thick consist mainly 
of halite. The upper surfaces of these zones might represent surfaces of desic- 
cation or near-desiccation that succeeded fairly long episodes in which saline 
lakes partly filled the basin. Thin layers of carbonates alternating with halite 
may represent shallow, short-lived lakes somewhat resembling playas. The 
“Mixed” layer as a whole implies warmer and drier conditions than those 
inferred from the Bottom Mud and Parting Mud. The climate represented seems 
to have ranged from one like that of today to a cooler, moister one, but with 
lower temperature not sustained as long as were the later conditions implied 
by the Bottom Mud and the Parting Mud. 

The events and climates inferred from the foregoing description of the 
stratigraphic units are rationalized in table 1. In the table the term nonpluvial 
refers to a climate similar to that of the present day, The term pluvial refers to 
a climate appreciably cooler and moister than that of today, No other signifi- 
cance is attached to these terms. 


rape | 
Rationalization of Stratigraphy and Shore Features of the Searles Basin 
with Climate. (Although the suggested correlation of stratigraphic units 
with strandlines seems probable, it is not established). 


Inferred 
Stratigraphic units r Implied events climate 
Basin floor 
Cessation of inflow; at 
Upper Salt ow strandline least near-desiccation; Nonpluvial 
possibly with interrup- 
tions 


Higher strandlines; Refilling of basin with 
Parting Mud upper limit unde- fresh water te an unde- Pluvial 
ermined termined altitude 


Cessation of inflow; at 

least near-desiccation; Nonpluvial 
possibly with interrup- 

tions 


Lower Salt Not identifie 


Higher strandlines, 
not separately identi- Filling of basin with 
Bottom Mud fied, Possibly partly fresh water to an unde- Pluvial 
or wholly obscured by termined altitude 
strandlines of the 
younger lake 


Basin at least partly 
filled with water much 
Strandline s, probably of the time, with occa: 
“Mixed” layer now obscured by sional shrinkage to 
strandlines of saturation point for Na 
later date carbonates, and possibly 
desiccation 


Nonpluvial to 
intermediate 
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RADIOCARBON DATES 
Dates from the Searles Lake sequence 

The discovery that the Parting Mud contains extractable organic matter 
opened the way for C'* dating of that layer, and made possible the identifica- 
tion and correlation, not only of the pluvial-lake phase represented by that 
layer. but also of the evaporite phases represented by the two salt bodies.° 
Samples of organic matter concentrated from drill cores were dated at the Uni- 
versity of Chicago in 1953 and 1954, and later others were dated at the labora- 
tories in Groningen and Yale Universities and in the U. S. Geological Survey, 
and at the Magnolia Petroleum Company Field Research Laboratory.® The 
samples dated are shown by number in table 2, which also shows the core from 
which each sample was taken and the substance used for dating. Figure 5 is a 
generalized column in which the dated samples are represented by black 
rectangles. The height of each rectangle shows the stratigraphic thickness for 
which the date is an average. The thicknesses vary with the concentration of 
C'* needed for processing in the various laboratories. From figures 3 and 4 it 
is seen that the thicknesses of the stratigraphic units vary slightly from one 
core to another within the central part of the Searles basin. Hence the column 
in figure 5, being generalized, includes a small distortion of the vertical dimen- 
sion, as samples from 7 different cores are plotted on it. With one exception 
noted hereafter, the positions of the samples shown in the figure are plotted 
according to their distance from the top or base of the stratigraphic unit in 
which they occur, whichever is the smaller. 

The rocks underlying the drainage basin that contributed to Ancient 
Searles lake include very few carbonate rocks. Hence it is likely that recycled 
ancient carbonates were not present in sufficient concentration to affect the C* 
dates by more than a slight amount. 

The most clearly established results of dating pertain to the Parting Mud, 
the unit on which dating efforts were concentrated first. The dates are con- 
sistent internally within statistics, and the results from the three laboratories 
represented are consistent also. As discussed in a foregoing section, the base of 
the Parting Mud shows that a dry or nearly dry saline flat in the Searles basin 
gave way to a lake in which mud was actively deposited. This implies that 
Owens Lake was then full to overflowing and was discharging actively into the 
Searles basin (fig. 1). Such an event demands a cooler and moister climate 
than that which preceded or followed the lake episode, as do also the pollen 
data for the Parting Mud, and some of the mollusks in the exposed silt on the 
assumption that silt and Parting Mud are correlative. 

According to the dates inflow of water began before 23,000 yr. B.P. (the 
average of the six finite dates from near the base of the Parting Mud), and 
ceased about 10,000 B.P. if cessation of mud deposition was caused by cessa- 
tion of overflow of Owens Lake. That inflow ceased because of climatic warm- 
° The dating possibilities were the subject of correspondence among W, A, Gale, Eliot 
Blackwelder, R. F. Heizer, Ernst Antevs, R. F. Flint, and W. F. Libby. 


* The samples marked C in table 2, column 1 were obtained from organic carbon extracted 
from pulverized core by means of organic solvents as described on a foregoing page, All 
the other organic-carbon samples listed were obtained from carbon recovered by burning 
pulverized core from which carbonates had been removed previously. 
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ing and drying is supported by the pollen statistics already mentioned, The 
Parting Mud therefore represents a pluvial episode at least 13,000 years long 
but not greatly in excess of that length. 

The rate of sedimentation can be approximated only crudely because each 
date is only the average date of a particular stratigraphic thickness of sediment. 
On the basis of the four finite Chicago dates (on organic carbon) the 5-foot 
thickness of mud measured from the center of sample C-897 to the center of 
C-896 was deposited in about 5,000 years, the next 3.5 feet, between C-896 
and C-895, in about 4.500 years, and the next 3 feet, between C-895 and C-894, 
in about 4,500 years. According to these rough values sedimentation rates were 
about 1 ft./1,000 yr. early in the period, decreasing to about 1 ft./1.300 yr. in 
the middle, and to 1] ft./1,500 yr. toward the end, averaging 1 ft./1.175 yr. for 
the whole period. 

On the other hand, the four Yale dates (Y-574b to Y-577b) on organic 
carbon’ give sedimentation rates of about 1 ft./500 yr. early in the period. 
about 1 ft./1,000 yr. in the middle. and almost 1 ft./3,000 yr. toward the end, 
averaging 1 ft./1.260 yr. for the whole period. It is not known which sequence 
of dates affords the more reliable basis of measurement; more detailed ex- 
amination of the sedimentary laminae might furnish a clue. 

From the dates of the Parting Mud it can be inferred that crystallization of 
the Upper Salt occurred during some part of the last 10,000 years. It was hoped 
that dates from the salt body itself would make possible greater precision, but 
the three values obtained have not accomplished this, Samples W-248 and 
RC-50c consisted of inorganic carbonate from trona. Their use was based on 
the principle that CO, dissolved in the lake, in approximate equilibrium with 
atmospheric CO., was the source for inorganic carbonates precipitated on the 
lake floor as well as for organic carbon. However, as figure 5 shows, not only 
are W-248 and RC-50c not consistent, but there is a small discrepancy, within 
statistics, between RC-50c and RC-50, organic carbon from the same 0.3 feet 
of sediment. Whether the latter discrepancy is real can be learned only from 
more dates, which are not available at present. A similar discrepancy appears 
in the group W-340, 341, and 343 from the Lower Salt, as noted hereafter, al- 
though the group GRO-1460, 1802. and 1805 agree within statistics, The dis- 
crepancies therefore may be related to contamination of the carbonates while 
in the ground or to laboratory procedures. Whatever their cause, for the present 
the two organic-carbon dates are preferred to the two carbonate-carbon dates 
from the Upper Salt. 

The strata beneath the Parting Mud are less well dated. The dates present- 
ly available indicate that the lake that earlier occupied Searles basin had al- 
ready existed for a long period at the time represented by the organic-carbon 
sample GRO-1464, about 46,000 yr. B.P. Mechanical sedimentation in the lake 
had ceased and precipitation of trona had begun prior to about 32.000 yr. B.P.. 
the approximate date of samples GRO-1805, 1460, and 1802 from the basal part 
of the Lower Salt. If we apply the mean rate of sedimentation indicated for the 
Parting Mud, about 1 ft./1,200 yr.. to the 10-foot thickness of sediment be- 


* These are considered more reliable than carbonate-carbon values, owing to the probably 
greater freedom of organic carbon from alteration by exchange after deposition. 
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Fig. 5. Generalized columnar section of upper part of sequence developed from cores 

central part of Searles basin, with positions and dates of C’* samples. Laboratory num- 

bers are keyed to table 1. Dates are quoted in years B.P. (before the present), Vertical 
scale in center column is not uniform with the two others. 


tween GRO-1464 and the top of the Bottom Mud, we get a theoretical elapsed 
time of 12.000 years. Subtracted from 46,000 years, the approximate date of 
GRO-1464, this gives 34,000 years as a tentative date of the top of the Bottom 
Mud. a value in fair agreement with 32,000 years for GRO-1460. Application 
of the average Parting Mud sedimentation rate is, however, based only on as- 
sumption; the rates for the two mud bodies could have been different. 

Although a date of about 32,000 B.P. is satisfactorily fixed for a zone near 
the base of the Lower Salt, samples from about the middle of that body give 
conflicting dates. Within a stratigraphic thickness of only 0.1 foot, organic 
carbon (W-343) gave ~29,500, inorganic carbon from carbonate (W-341) 
gave ~ 23.000, and wood from a fossil twig (W-340) gave ~26,700 years. The 
inorganic-carbon date is rejec ted as being inconsistent with the dates of GRO- 
1802, 1805, and 1460 and patently too young. The conflict between RC-50 and 
RC-50c, in the Upper Salt body, is analogous. The other two (W-340, W-343) 
agree within statistics and their average, ~ 28,100, is accepted tentatively as 
the best value available. 

This value is not necessarily inconsistent with the adjacent GRO dates 
shown in figure 5, because two cores, 3.75 miles apart, are involved, In one of 
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them (GS-27), located near the margin of the basin, the Lower Salt is only 
about 23 feet thick. whereas in the other (X-20), near the center of the basin. 
the Lower Salt is 39 feet thick. Therefore, although W-343 and W-340 are at 
almost the same position measured from the base of the Lower Salt as is GRO- 
1802, they are probably younger than GRO-1802. On a basis of the thicknesses 
of the Lower Salt in the two cores, the position of the W samples has been 
projec ted into the generalized column for the central part of the basin as indi- 
cated in figure 5. Although determined mechanically this position is believed 
to be more realistic than one based on measurement from the base of the Lowe 
Salt. Crystallization of the Lower Salt must have been embraced within the 
interval >32,000 B.P. to somewhat more than 23,000 B.P., a period probably 
more than 9,000 yr. long, and desiccation or near-desiccation of the Searles 
basin at the end of Lower Salt time then occurred shortly before 23,000 yr. 
B.P. The duration of that condition seems to have been short. 

We can attempt to estimate crystallization times for the Lower and Upper 
Salt bodies in another way. The thickness of the Lower Salt is a crude measure 
of the amount of dissolved substances brought into the earlier lake after it 
ceased to overflow, and is a still cruder measure of the time elapsed between 
cessation of overflow and near-desiccation, if we assume that the annual incre- 
ment is known and was constant. and that the rate of evaporation also was 
constant. The C’* dates imply that about 13,000 years were required to bring 
enough dissolved substances into the Searles basin to precipitate, afterward. 
the (average) 80-foot layer (2.4 billion tons) of salts constituting the Upper 
Salt. This is at a rate of 1 foot (30,000,000 tons) of precipitated salts per 162 
years of antecedent contribution. Assuming the areas of the Upper Salt and 
Lower Salt are equal and assuming a similar rate of contribution, the (average } 
30-foot (0.9 billion tons) Lower Salt would have required 4,860 years of ante- 
cedent contribution after any overflow from the Searles basin had ceased. 

Similarly, evaporation of Searles Lake to dryness during precipitation of 
the 80-foot Upper Salt (possibly with interruptions by freshening) required 
a maximum of 10,000 years. At an assumed similar rate of evaporation. and 
neglecting interruptions, precipitation of the 30-foot Lower Salt would have 
required a maximum of 3,750 years. Obviously these calculations are subject 
to large errors because of the assumptions made, and because they neglect loss 
of solutes through possible overflow from the Searles basin. In any case they 
refer only to the accumulation of the evaporites; they do not take account of 
the time interval between completion of precipitation of the Lower Salt and 
restoration of the lake represented by the Parting Mud. For what they are 
worth, the values are compatible with the C™ dates. 

We have, then, a record, controlled by C** dates, of two pluvial climates 
marked by deep-lake phases, each followed by a nonpluvial climate. The later. 
existing nonpluvial climate is marked by desiccation; the earlier one was 
marked by a similar condition or by something approaching it. The earlier 
pluvial began long prior to 46,000 B.P. and was on the wane before 32.000 
B.P.; the later pluvial began at or before 23,000 B.P. and was on the wane 
before 10,000 B.P. The earlier time of maximum dryness immediately preceded 
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the onset of the later pluvial. The later desiccation is still in effect today; how 
long ago it started is not known. 


Date of last overflow of Owens Lake 

H. S. Gale (1915, p. 260-264) calculated the “age” of Owens Lake by 
dividing the amounts of chlorine and sodium in solution in the lake in 1912 
by the respective annual inflow of these elements estimated from sampling of 
Owens River. By “age” he meant the time since Owens Lake ceased to over- 
flow and began to accumulate these elements. He arrived at values of 4,200 
years (C1) and 3,500 years (Na).* 

However, Antevs (1938, p. 191) showed that these values could be ex- 
plained alternatively as representing the time since a rebirth of Owens Lake 
following evaporation to dryness, which should have occurred in the period 
(Hypsithermal and its equiv alénts) centered around 5,000 yr. B.P. 

Whether H. S. Gale’s values do in fact measure the cessation of Owens 
Lake overflow is not clearly evident from the stratigraphy of the Searles sedi- 
ments. Those sediments tell us with certainty only that the approximate date 
10,000 yr. B.P. sets a maximum limit to the time when the lake water became 
saturated with trona. It can be argued that that date also approximately meas- 
ures the time when long-continued inflow to Ancient Searles lake, and hence 
outflow from Owens Lake, ceased. Other sources of inflow to Ancient Searles 
lake are negligible; water to maintain the lake would have had to come from 
a much greater segment of the eastern watershed of the Sierra Nevada than 
that which now furnishes water to Indian Wells Valley. However, the Upper 
Salt contains several mud partings. Conceivably these represent re-freshenings 
of Searles Lake, brought about by overflows of short duration from Owens Lake. 
Such overflows could have transferred accumulations of salts in solution from 
Owens Lake to Searles Lake. It is possible therefore that temporary overflows 
from Owens Lake to Searles Lake might have occurred less than 10,000 years 
ago. On the other hand if the mud partings represent sediment washed into the 
center of the basin from local sources, we know of no evidence that Owens Lake 
overflowed after the Parting Mud ceased to be deposited. 


Possible dating of sediments in other basins 

Probably sediments of other lakes of the ancient chain likewise contain 
extractable organic matter. If so, and if suitable samples could be secured, the 
sequence of events inferred for Searles Lake could be checked in part by 
stratigraphic study and C'* dating of such substances in Death, Panamint, and 
Owens Valleys. Borings made in Death Valley have established the presence, 
beneath its floor, of alternating layers of clay and evaporites (H. S. Gale, 
1914b, p. 412-413). Dates on clay layers could be compared with the inference 
drawn by Blackwelder (1933, p. 470) from geomorphic features, that Lake 
Manly existed in Tahoe time but not in Tioga time. Borings made in Panamint 
Valley (Smith and Pratt, 1957, p. 51-62) show clay and silt with minor gyp- 
sum, carbonates, and halite. If dates are possible they should throw light on 
the time of last overflow of Ancient Searles fake. 


An early printing of H. S, Gale’s report, still existing in some libraries, gives much larger 
values, which are corrected in a later printing. 


‘ 
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CORRELATION OF EVENTS IN THE SEARLES BASIN WITH EVENTS ELSEWHERE 
Glacial events 

Inclusive dates of the later pluvial period inferred at Searles Lake 
( >23,000 to ~ 10,000 B.P.) approximate the inclusive dates of the classical 
Wisconsin glaciation of central United States. In that glaciation the Laurentide 
Ice Sheet invaded the Great Lakes region about 25,000 B.P., reached a maxi- 
mum at or later than 18.000 B.P.. and after lesser readvances began its final 
retreat about 10.700 B.P. (See Flint and Rubin, 1955; Flint, 1955, 1956). 
The agreement of the two pairs of inclusive dates indicates a nearly 1-to-1 time 
relationship between the classical Wisconsin glaciation in the eastern part of 
the United States and the pluvial climate represented by the Parting Mud. This 
time coincidence agrees with stratigraphic evidence of coincidence of glacial 
with pluvial events elsewhere in the Basin-and-Range region (Gilbert, 1890, 
p. 305-311: Hunt, Varnes, and Thomas, 1953, p. 28; Putnam, 1950), and 
thereby indicates that the last glacial and the last pluvial, at least within the 
belt of westerly winds on the North American continent, are sensitively related 
to a climatic cause common to both. 

The close time agreement between the Parting-Mud pluvial and the classi- 
cal Wisconsin glaciation implies that the Bottom-Mud pluvial should find a 
correlative in the glacial sequence in North America. Evidence of glaciation in 
the range 40,000-50,000 yr. B.P. should be looked for in that sequence. 

There exists no physical connection on the basis of which sediments in the 
Searles basin can be correlated with glacial features in the Sierra Nevada. 
However, the similarity of their climatic implications and stratigraphic posi- 
tions suggest tentative correlation of the Bottom Mud with the Tahoe glacial 
unit and of the Parting Mud with the Tioga glacial unit. The two glacial units 
(Blackwelder, 1931) are not far apart in time and represent the two last major 
slacial climates in the Sierra Nevada, The correlation, which seems reasonable, 
could be tested best hy C'* dating. 


Lakes Bonneville and Lahontan® 


The dated sequence beneath Searles Lake can be considered profitably in 
connection with the record of Ancient lakes Bonneville and Lahontan. Since 
the classic work of Gilbert (1890) on Lake Bonneville, geologists have thought 
of the high-water phases of that lake as approximating glacial maxima, and of 
the low-water or desiccation phases as recording glacial minima; the Searles 
Lake data support that opinion, as do the minor recent fluctuations of Great 
Salt Lake, which are correlated with climatic variations, The Searles data like- 
wise are compatible with the dates on core material from between 18 and 30 
feet below the bottom of Great Salt Lake, and consisting of sediments of the 
last freshwater phase of the lake. The dates are 16,850+300 (W-321, organic 
carbon) and 16,680+300 (W-319, inorganic carbonate) (Rubin and 
Alexander, 1958, p. 1484). 

The Searles Lake dates suggest further that net evaporation was in prog- 
ress during a substantial later portion of the Parting Mud phase, because by 
the end of that phase, around 10,000 yr. B.P., the concentration of sodium 


* A pertinent and informative paper (Broecker and Orr, 1958), not published until after 


the completion of the present paper, is recommended as a reference. 
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carbonate in the lake water had reached the precipitation point for trona. These 
dates are compatible with a group of C'* dates (discussed by Jennings, 1957) 
from Danger Cave, Utah, in the western part of the Bonneville basin only a 
little more than 100 feet above Great Salt Lake. The dates of 6 samples from 
1 culture layers in the cave are internally consistent. The oldest, approximately 
11.500 yr. B.P., directly overlies the beach gravel of a late phase of Lake 
Bonneville. The stratigraphy shows that the site was not submerged after the 
layer containing the sample was deposited. Hence Lake Bonneville was under- 
going net evaporation, and uncovered the cave, at some time prior to 11,500 
B.P.. and after that date did not rise to the altitude of the cave. The compati- 
bility of these Bonneville and Searles dates suggests that variations in the two 
lakes may have been broadly synchronous. 

The relations at the Leonard Rockshelte’, Nevada (Heizer, 1951) re- 
semble those at Danger Cave. Situated in the downstream part of the Hum- 
holdt Valley, the rockshelter contains Lake Lahontan beach gravel at elevation 
1.175 feet, more than 150 feet above the floor of Humboldt Dry Lake. Although 
not specified, it appears from the facts published that after the beach gravel 
was deposited the rockshelter was not later submerged. The date of bat guano 
immediately overlying the beach gravel is about 11,200 years, very close to 
that of the lowest culture layer at Danger Cave. 

Compatible with those dates is the date of wood from the lowest culture 
layer in Fishbone Cave, facing Winnemucca dry lake, also in the Lake 
Lahontan basin and at an altitude less than 100 feet below that of the dated 
layer at Danger Cave. The site is reported to have been occupied by man 
shortly after the cave emerged; the date is approximately 11,200 years 
(Broecker, Kulp, and Tucek, 1958, p. 163). 

According to the C'* dates the beach at Leonard Rockshelter and that at 
Danger Cave, formed during fluctuating subsidence of the Bonneville and 
Lahontan pluvial lakes, were built at some time prior to that of the Two Creeks 
peat in eastern Wisconsin. That peat was overwhelmed by glaciation, the 
Valders readvance, believed to have reached its maximum around 10,700 B.P. 
(Flint, 1956). If the climatic change that caused that readvance was reflected 
in the regimens of Lakes Bonneville and Lahontan, any resulting rise of the 
lake surfaces was limited, for it failed to reach the altitudes of Leonard Rock- 
shelter and Danger Cave. The change has not been discerned in the stratigraphy 
of Searles Lake. Apparently therefore it was not conspicuous in the southern 
part of the Basin-and-Range region. 

If glacial readvances such as the Valders were responses to local variations 
in a general continent-wide rise of atmospheric temperature, they are expect- 
able where a large ice sheet is present in a state of metastability, But it is 
doubtful whether we should expect major response in the Basin-and-Range 
region, where glaciers were small and where, with the waning of Wisconsin 
glacial climate, the atmospheric circulation pattern may already have shifted 
enough to remove that region from the belt of “pluviality.” 

The fact that the sandals found at Fort Rock Cave, Oregon, and dated at 
9.000 B.P. (ref. in Flint and Deevey, 1951, p. 278). were woven from sage- 
brush is compatible with the view that shortly after the time of the Valders 
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maximum even the extreme northern part of the Basin-and-Range region was 
becoming arid. 

Although the Searles Lake dates and inferred events can be reasonably 
correlated also with dated events recorded in parts of Europe and beneath the 
deep-sea floor, in the present paper no attempt is made to suggest such wider 
correlation. 


SUMMARY 

The Searles basin was one of a chain of four basins that held temporary 
Pleistocene lakes. It received influent drainage only when the Owens basin 
overflowed; this could have occurred only when the climate was cooler and 
moister than today’s. Such a former climate is supported by fossil aquatic 
mollusks and by fossil pollen. 

The presence of numerous strandlines around the basin establishes that a 
lake formerly existed, but it has not been possible to correlate the strandlines 
with the lake-floor sediments known from core borings. 

At least two former lakes of different ages are indicated by two pairs of 
strata revealed in the cores. Each pair consists of clay and silt overlain by 
evaporites. The earlier lake lasted from before 46,000 to nearly 23,000 B.P., 
by which time it was dry or nearly so. The later lake formed about 23,000 B.P. 
and ceased to receive inflow about 10,000 B.P. The indicated fluctuations are 
believed to have been related directly to fluctuation of the regional climate. 

Data from the Searles basin do not yet afford a decisive basis for judg- 
ment on the validity of H. S. Gale’s calculation of the “age” of Owens Lake, 
nor on the times when Ancient Searles lake overflowed. The history and 
chronology of the later lake are compatible with published data on the history 
of Lakes Bonneville and Lahontan. 


The inclusive dates of the later lake in the Searles basin approximate 


those of the classical Wisconsin glaciation. Hence they confirm the synchrony 
of the last glacial with the last pluvial climate within the belt of westerly winds 
over the North American continent. 
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GLACIAL FLUTINGS IN CENTRAL 
AND NORTHERN ALBERTA! 
C. P. GRAVENOR and W. A. MENELEY 


ABSTRACT. An examination of vertical air photographs covering about 175,000 square 
miles of northern and central Alberta has revealed the presence of large fields of glacial 
lutings. The flutings show two dominant trends, one in a south to southwest direction and 
the other in a southeast direction. Although the majority of the flutings are formed in till, 
n northeastern Alberta the flutings are erosional in character and developed in Precam- 
brian granite gneiss. Striae on the Precambrian flutings are roughly parallel to the ridges. 
ill fabric determinations made on the southeast flutings of the North Battleford, Sas- 
katchewan district show that there is a marked change in fabric with depth in the ridges. 
ng development apparently bears little relationship to topography, It is believed that 
flutings originated from alternating high and low pressure zones at the base of the ice. 


INTRODUCTION 

The word “fluting” has been used to describe parallel shallow grooves on 
till plains (e.g. Chapman and Putnam, 1951, p. 15). In this report the term 
tluting is used in its widest possible sense and except where specified includes 
low-lying drumlins and giant grooves as well as parallel shallow grooves. The 
fact that drumlins—of a variety of form—and flutings are gradational (plate 
1) has been pointed out by several geologists (e.g. Smith, 1948; Gravenor, 
1956; and Lemke, 1958). 

An examination of vertical air photographs and mosaics of Western 
Canada shows that flutings are very common and in the absence of striae and 
well-developed end moraines, provide a good clue to the direction of glacier 
movement. There are two major ice-movement directions in northern and 
central Alberta, one to the south and southwest and the other to the southeast. 
The significance of these two directions is as yet rather obscure but there is 
preliminary evidence to suggest that certain flutings which trend in a south- 
easterly direction were created by minor readvances of the retreating ice. 

Most of the flutings examined are composed of sandy to clayey till except 
n northeastern Alberta where they are formed in Precambrian granite gneiss 
and porphyritic granite. Unfortunately exposures in fluting are extremely rare 
and the possibility exists that many of them contain incorporated stratified 
material. 


FLUTING DIRECTIONS IN NORTHERN AND CENTRAL ALBERTA 

Vertical air photographs covering about 175,000 square miles of northern 

and central Alberta have been examined and fluting directions recorded, The 
results of this study are shown on figure 1, Other glacial features, such as 
crevasse fillings and a variety of minor end moraines, can be confused with 
poorly developed flutings and consequently such doubtful features have been 
carefully excluded from this study. Undoubtedly as more detailed mapping of 
surficial deposits is carried out in northern and central Alberta more ice-move- 
ment directions will be added. Most of the measurements were made on the 
shallow-groove type of fluting but in some instances—for example in the area 
Paper presented before the Canadian Association of Geographers, Edmonton, Alberta, 
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Most of the ice-movement directions shown on figure 1 were developed 
from continental glaciers which advanced in a west to southwest direction from 
the Keewatin ice center. An exceplion is found in the area west of Edmonton 
between the North Saskatchewan and Athabasca Rivers where it is believed 
that the flutings resulted from Cordilleran ice which flowed in a northeast 
direction out of the foothills and then swung to the southeast more or less 
parallel to the foothills of the Rocky Mountains. 

The majority of the ice-movement directions over northern and central 
Alberta are towards the south and southwest and are in reasonably good 
agreement with the ice-movement directions predicted from the use of indi- 
cators (Gravenor and Bayrock, 1956). There are, however, a significant num- 
ber of directions toward the southeast (Plate 1B)—particularly in the Lac la 


Biche-Vermilion area. The Lac la Biche flutings extend in a southeasterly di- 
rection to a point 30 miles southeast of North Battleford, Saskatchewan, From 
the fact that southwest directions exist north and south of the Lac la Biche 


southeast trend it might be suspect ted that these two sets of directions are of 
different age. There is some data—which is presented under fabric studies 

which suggests that the Lac la Biche southeast flutings were probably developed 
during deglaciation of east-central Alberta and have no relationship to the 


direction of regional ice advance. 


LOCAL DISTRIBUTION OF FLUTINGS 

Effect of topography.—Although topography may have been an influence 
on the direction of ice movement in central and northern Alberta the relation- 
ship of topography to fluting development is rather obscure. In the Birch 
Mountains (plate 2B) 70 miles northwest of McMurray (Ells, 1932) and in 
the Naylor Hills 100 miles northwest of Peace River the flutings are developed 
only on the tops of highlands. In the Red Deer district (plate 2A), on the other 
hand, the flutings are best developed in a sidehill position of a preglacial 
valley. 

In the Stettler district (fig. 2) the flutings were developed from the glaciet 
ice which moved uphill in a southwesterly direction. The fan-like distribution 
of these flutings suggests that the ice spread laterally as it was advancing to the 
south. It is difficult to visualize how such a distribution of flutings could origin- 
ate from a regional ice flow and it is suggested that the Stettler flutings were 
created by a minor readvance of the glacier. The last glacier probably retreated 
well to the north of Stettler and then readvanced, locally, sending a small ice 
lobe into the Stettler district. It is interesting to note that linear trends visible 
in the hummocky dead-ice moraine five miles north of Stettler show the same 
orientation as the flutings. Although the age of the morainal ridges is post- 
fluting it appears that the same stress distribution in the ice which gave rise 
to the flutings in part controlled the positioning of some of the dead-ice fea- 
tures, 

In summary, flutings are found on the tops of highlands, on sidehill posi- 
tions in lowlands and on uphill positions, Thus it appears that while the ice 
movement may have been affected by topography, fluting development is in- 
dependent of topographic control. A similar conclusion has been reached by 
Smith (1948) in his study of giant glacial grooves. 
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EROSIONAL FLUTINGS ON THE PRECAMBRIAN 
The area north of Lake Athabasca in the extreme northeast corner of 
Alberta is underlain by massive granites and metasediments of Precambrian 
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ige. The surface of this portion of the Canadian Shield is gently rolling and 
broken, locally, by ice-scoured valleys and fault scarps. In the Andrew Lake 
district (fig. 3) glacial erosion on the Precambrian rocks produced a field of 
flutings which trend north 60 degrees east. To the east in northwestern 
Saskatchewan the erosional flutings are gradational to a drumlin and fluting 
field. The flutings in the Andrew Lake area transect the bedrock structures 

hg. 3) and hence differential erosion through variation in lithology has been 
significant. In many instances feldspathic pegmatite dykes which cut across 
-everal flutings have been eroded to the same depth as adjacent softer meta- 
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These erosional flutings are similar in size and shape to the giant grooves 
developed on sedimentary rocks in the Mackenzie district—described by Smith 
(1948). In general, the fluting ridges lie about 10 to 25 feet above the adjacent 
grooves and the individual ridges are spaced about 300 feet apart. Some of the 
fluting ridges can be traced for a distance of several miles, but most are con- 
siderably shorter. 

Many hundreds of stria measurements made on the sides and tops of 
Andrew Lake flutings show a parallelism to the ridge orientation (fig. 3). Thus 
it is reasonably certain that the flutings were developed parallel to the direc- 
tion of ice movement, although it should be kept in mind that the striae were 
developed at a later date than the flutings and perhaps under different ice 
conditions to those which actually created the flutings. The fact that the 
Andrew Lake flutings are gradational to drumlins and are similar in size. 
shape and spacing to till flutings would suggest that the erosional process is an 
important factor in the development of flutings of all types 


FLUTING WAVELENGTH 

The term fluting wavelength as used in this report is defined as the dis- 
tance between the crests of adjacent ridges or troughs of adjacent grooves. 
measured perpendicular to the long direction of the flutings. Fluting wave- 
length measurements were made on air photographs along lines drawn two- 
thirds to one mile apart normal to the fluting trend. In most cases wavelength 
measurements were made between adjacent ridges but in some cases the 
grooves were used in order to get the best values. The wavelengths of all visible 
flutings which crossed each line of measurement were recorded. Wavelengths 
greater than one mile are rare and are not reported. 

Wavelengths of flutings in five districts of Alberta and Saskatchewan were 
measured. The results were plotted as histograms (fig. 4) with the wavelength 
the abscissa and the frequency of occurrence the ordinate, Statistical tests show 
that the histograms have a significant deviation of wavelengths from a normal 
distribution. There is a preferred wavelength of 300 to 400 feet present in all 
cases and a secondary preferred wavelength of 600 to 700 feet in all cases 
except for the North Battleford flutings. The secondary peak is best developed 
on the Cold Lake flutings and may be related to the fact that a portion of these 
flutings is composed of elongate drumlins (plate 1B). 

The similarity of fluting wavelength developed under widely differing 
conditions of topography and lithology suggests that wavelength is essentially 
independent of topographic or lithologic control. The wavelength must there- 
fore be controlled by some physical property of the glacier ice which gives rise 
to periodic variations in erosive capacity in a direction transverse to the direc- 
tion of flow. Until more is known about the physics of ice flow it is impossibl 
to suggest a reason for periodic variation of the erosive capacity in ice. 


MICROFABRICS~ 


Macrofabric studies made on drumlins generally show a prominent linea- 
tion parallel to the long axis of the drumlin (e.g. Donner and West. 1955: 


* The method of preparation, analysis and interpretation will be the subject of a further 
paper which is now in the preparation stage. 
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Wright. 1957). Due to the absence of deep exposures it was decided to study 
the fabric of flutings by microfabric techniques. According to Seifert (1954) 
and Harrison (1957) microfabric studies give similar results to those obtained 
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by macrofabrics and hence appear to be a reliable technique in the study of 
till fabric. 

Microfabric analyses were carried out on oriented undisturbed samples o! 
till obtained from auger holes in a fluting at North Battleford, Saskatchewan. 
Shelby spoon samples of three-inch diameter were used to take samples at in 


} 


tervals below the crest of a fluting ridge. The samples were impregnated wit! 
a thermosetting plastic and a minimum of three mutually perpendicular thin 
sections was cut from each sample. The least projection elongation direction 
(Dapples and Rominger, 1945) of the apparent shape of about 200 particles 


e bservat 


Fig. 5. Microfabric diagrams (radius of circle = 10 observations). 
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was measured in each thin section and plotted on rose diagrams, Microfabric 
diagrams of two typical samples are shown in figure 5. 

In order to define the attitude of any particle in space it is necessary to 
measure the orientation of two shape parameters—the long axis and the plane 
of maximum projection area. In this study the preferred orientation of the long 
axis and the plane of maximum projection area are described by the standard 
structural symbols for lineation and foliation respectively. 

The interpretation of the two-dimensional orientation data is based upon 
graphic analysis of the relationships between the two-dimensional apparent 
shape in thin section and the long axis and maximum projection area of the 
particles. If the foliation dips less than 45 degrees then the preferred orienta- 
tion determined from a thin section cut in a horizontal plane approximates the 
strike of the lineation. The preferred orientation directions determined in two 
vertical thin sections may be considered as the apparent dips of the plane of 
foliation. The strike and dip of the foliation and the plunge of the lineation 
may then be determined by the methods described by Billings (1954, p. 442- 
143). 

Ata depth of four to five feet below the crest of the fluting ridge, micro- 
fabric determinations (fig. 5) show a prominent lineation which strikes parallel 
to the long direction of the fluting—northwest—and plunges northwest at 20 
degrees. The foliation dips west at 30 degrees. A similar microfabric is present 
at a depth of 10 feet with the exception that lineation plunges at a very small 
angle—two degrees—to the northwest and the foliation dips eight degrees 
southwest. 

At depths greater than 10 feet the microfabric changes markedly although 
there is no apparent change in the mechanical composition of the till or no 
evidence to suggest that more than a single till is present. At 11 feet the micro- 
fabric (fig. 5) has a lineation which strikes a few degrees east of north and 
plunges at 25 degrees to the north. The foliation dips north-northwest at 27 
degrees. At 13 feet the microfabric (not figured) has only a lineation striking 
north, 

The till constituting the upper 10 feet of the fluting ridge is interpreted 
as till which was eroded from the neighboring grooves. The till was transported 
obliquely from the grooves and deposited to form the ridge by ice advancing 
parallel to the trend of the flutings. It is believed that the till was not trans- 
ported far before being deposited on the ridges. 

The North Battleford flutings on which this microfabric study was made 
form a continuation of the southeasterly-trending flutings of the Lac la Biche 
area of Alberta (fig. 1). It is interesting to note that the south to southwest 
lineation found below a depth of 10 feet in the North Battleford flutings is 
almost parallel to the southwest-trending flutings found north (Cold Lake) 
and south of the North Battleford-Lac la Biche flutings (fig. 1). Although more 
data are required it might be speculated that the southwest Cold Lake flutings 
were developed by a regional ice advance* and later during deglaciation the 
southeasterly Lac la Biche trend was produced by a local readvance or succes- 


* This hypothesis is supported by the regional southwest fluting and drumlin trends over 
northern Saskatchewan shown on the Glacial Map of Canada (1958). 
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sive readvances of the ice—perhaps flowing parallel to the retreating ice front. 
If such is the case the data cast some doubt on the validity of flutings as in- 
dicators of regional ice-movement direction and suggest that fabric determina- 
tions made at depth might prove more reliable for this purpose. 


ORIGIN OF FLUTINGS 

It is not the purpose here to discuss the many hypotheses which have been 
proposed for flutings and allied features. As more data are collected, however, 
it becomes apparent that erosion has played a significant role in the develop- 
ment of such forms. The fact that flutings are well developed in sedimentary 
and metamorphic rocks and that such flutings are gradational to drumlins and 
flutings developed in bedrock or in sandy till make it certain that flutings and 
drumlins can be produced through erosion. 

In dealing with such materials as clayey till, factors other than erosion 
may have a bearing on the production of flutings. Clayey till, when water 
saturated, is highly plastic and can be moulded and shaped quite readily. If 
unfrozen, water-saturated clayey till exists at the base of a glacier it could 
readily be pressed into tunnels or other openings at the base of the ice. Dyson 
(1952) has suggested that small-scale flutings could be formed in this way. 
The regularity of spacing, the parallelism and length of the flutings found in 
Western Canada all make it appear unlikely that tunnels or other basal-ice 
openings were responsible for fluting development. At least it is difficult to 
visualize how such openings could be created at regular intervals and main- 
tained in an oper state under a thick, moving ice cover. 

Lemke (1958) believes, with some reservations, that those flutings (nar- 
row linear drumlins) which contain stratified materials, “ . .. were formed by 
glacier ice advancing over stratified deposits, eroding and shaping the deposits 
into ridges and then plastering a layer of till over the deposits”. With the ex- 
ception of the layer of plastered till the explanation given by Lemke is quite 
similar to that suggested by Gravenor (1953) for the origin of drumlins, 


Any proposed origin must take cognizance of the following factors. 

(1) Flutings and drumlins are gradational features 

(2) Flutings are formed of till, stratified drift and bedrock 

(3) Flutings are regularly spaced, show well-developed parallelism and 
are quite long (up to several miles) 

(4) There is some evidence to suggest that many flutings were created 
by a local readvance of a glacier which overrode existing glacial 
deposits 

(5) Preliminary fabric data show that there is lineation parallel to the 
ridges and foliation which dips away from the ridge tops. 


From the fact that flutings are found in bedrock and in some cases con- 
tain stratified materials (Lemke, 1958) the writers believe that erosion has 
played a significant role in the development of such forms. The preliminary 
fabric data suggests that till which is found on the tops and sides of flutings 
was placed there, partly, through the component of ice flow. perpendicular to 
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the long trend of the ridges. Considering both of these factors—erosion and 
fabrie—a picture can be built up of an actual process of formation. It is sug- 
gested that flutings have been developed through the action of alternating 
parallel high and low pressure zones at the base of the ice. Figure 6 shows a 
hypothetical cross-section through the base of a glacier perpendicular to the 
direction of ice flow. Material eroded from below the high pressure zones is 
incorporated with the ice and moves downstream in a curved path toward the 


ICE 


Fig. 6. Hypothetical cross-section through the base of a glacier to illustrate a pos- 
sible mode of formation of flutings, Arrows indicate the lateral component of glacier flow. 


low pressure zone on the sides and top of the fluting ridge. By this mechanism 
the ice between the grooves becomes debris-clogged and till is deposited over 
the core of undisturbed material. This mechanism presupposes that flutings 
were created from pre-existing glacial deposits—an assumption which must be 
made in any erosional hypothesis and is proven in the case of flutings formed 
in bedrock. 

It is believed that this mechanism of fluting formation acounts for the 
erosional character of flutings, diverse composition and observed fabric. It 
might be argued that instead of the till being eroded from the grooves it was 
simply squeezed from the high pressure to the low pressure zone, Such a 
process necessitates unfrozen water-saturated till at the base of the ice, While 
this is a possibility it nevertheless enters another assumption which is un- 
warranted at this time. Additional work on microfabrics is now in progress at 
the Research Council of Alberta and it is hoped that such studies will result 
in a clearer picture of fluting development. 
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A FURTHER ATTEMPT TO DETERMINE 
THE VAPOR PRESSURE OF BRUCITE 


W. S. FYFE 


ABSTRACT. An attempt has been made to determine the equilibrium vapor pressure 
of brucite by studying pressure-volume relations at constant temperature. The results ob- 
tained agree closely with those found previously by Kennedy but differ from those found 
by Roy and Roy. No reasonable explanation of this difference is offered except by the 
general observation that quenching methods are fundamentally sound only when the quench 


preserves the phases present under the conditions of the experiment. 


Recently in this journal Roy and Roy (1957) have presented data on the 
vapor pressure of brucite which confirms the earlier results found by Roy, 
Roy and Osborn (1953). Their data obtained by a synthetic method conflicts 
with the more direct determination made by Kennedy (1956). Kennedy’s data 
are very close to the results expected from available thermodynamic data as 
evaluated by MacDonald (1955). 

Some measurements have recently been carried out in this laboratory 
and, in view of the discrepancies in other results, are reported here. The 
method used was, in principle, similar to that used by Kennedy and depends on 
a study of pressure-volume relations in the system. 

Anhydrous MgO (5-12 gms.) was packed into a pressure vessel of about 
11 ce internal volume. The vessel was heated to temperature in an electric 
furnace and thermocouples were placed in holes in the front, rear and center 
of the vessel. When the temperature was steady, water was slowly introduced 
into the vessel using a compressor of a type similar to that used by Kennedy 
(1950) in P-V-T measurements. The volume of water introduced, or the piston 
displacement, could be accurately measured. Pressures were measured with 
Bourdon gauges and thus pressure-volume curves could be constructed. 

lemperature readings on the three thermocouples indicated that control 
and variation was within + 3°C. The absolute accuracy of pressure measure- 
ments is probably no better than + 25 bars although the reproduceability of 
readings is within + 7 bars. The gauges used checked well against a manganin 
coil gauge and there is no reason to believe that there are any serious errors. 

If the system behaved ideally we would expect a pressure-piston displace- 
ment curve, both for increasing and decreasing displacements, of the type 
shown in figure 1. The plateau on this curve would represent the equilibrium 
pressure at the particular temperature chosen. In an actual experiment this 
ideal form was not observed due to the rather slow kinetics of the reaction in 
the region near the equilibrium pressure. The technique adopted was to intro- 
duce constant increments of water and then follow the pressure variation. Thus 
the actual results appeared as in figure 2 (curve A). If there were serious 
gradients in the sample the plateau region would be curved but no serious 
curvature was observed. 

After hydration was complete the piston was withdrawn in stages and the 
behavior observed was as in curve B (fig. 2). A total cycle normally occupied 
about seven days as times to reach equilibrium are quite large. 

It will be noticed that the equilibrium values recorded for the hydration 


729 


730 W. S. Fyfe 


step are above those for the dehydration reaction. This gap (fig. 4) was 
normally of the order of 100 bars. If all the reactions involved in the process 
were fast there should be no such gap and some more detailed work was car- 
ried out on the rates. It was hoped that a rate equation might allow extrapola- 
tion to infinite time but this is not possible. After 12 hours the pressure changed 


A A 
AA 
| 
| | | 
/ 
Fig. 1. Curve illustrating Fig. 2. Actual behavior observed 
ideal behavior during hy- during hydration of MgO (curve 
dration of MgO at constant A) and dehydration of Mg(OH) 
temperature and volume. D (curve B). 


is the piston displacement. 


extremely slowly, around 7 bars in 24 hours. Such values could not be read 
with sufficient accuracy to allow any suitable extrapolation. Although the de- 
tailed shape of the short time portion of these curves was not very reproduce- 
able the 24 hour values were constant to within + 7 bars. Two typical rate 
curves for hydration and dehydration are shown in figure 3. 

On account of this difficulty in the results summarized in figure 4 curves 
are drawn for points from hydration (A) and dehydration (B) and it is as- 


2 4 6 8 10 12 14 16 18 20 
Time hours. 


Fig. 3. Curves illustrating the rate of hydration of MgO (A) and dehydration of 
Mg(OH). (B) at 542°C. 
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sumed that the true curve (C) is intermediate between the two. The points are 
drawn to indicate their uncertainty. 


T 

Bars. 

800 - t 

600 - + 

500 |; + 

400 + t 

200 
500 520 540 560 580 


Fig. 4. Measured vapor pressures for hydration (A) and dehydration (B). The 
equilibrium curve is assumed to lie between (C), Curve D is that obtained by Kennedy 
(1956). 


DISCUSSION 

There is no need for elaborate discussion of this paper. Curve D of figure 
b represents Kennedy's values, and it will be noticed that the results obtained 
here are very similar. It is difficult to see any mechanism which could lead to 
these results being in error by the amount required to match the data of Roy 
and Roy (1957). The following comments on their work are perhaps pertinent: 
(a) It appears that water was present during heating and equilibrium was 
approached from one direction only. 
(b) Their point I (p. 577) states: “If in a p-t quenched run of the usual 
length all the samples yielded well-crystallized brucite the point MUST be in 
the brucite field”. This is a true statement but no results are recorded in the 
table which fulfil th condition. Basically, the phases present during the run 
afe recognized not by the chemical composition or structure of the products, 
but by microscopic texture. This requires skill and experience and there is no 
doubt that the results may be correctly analysed. De Boer (1957) has discussed 
orientation of product lattice on parent lattice and while this is true for de- 
composition it may also be true for formation and make interpretation difficult. 

The difference between the hydration and dehydration steps found in the 
present work requires some explanation, but the writer feels that surface effects 
caused by grain size could be adequate. The observed differences could be 
caused by a free energy change of 200 cals. and this is a small term compared 
with the extreme differences discussed by MacDonald (1955). 
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NEUTRALIZATION CURVES AND THE FORMULATION 
OF MONOVALENT CATION EXCHANGE PROPERTIES 
OF CLAY MINERALS* 


PALL D. BLACKMON 


ABSTRACT. Exchange constants derived from observed titration curves for three clay 
minerals (beidellite, “illite”, and bentonite) were used to calculate new titration curves 
against potassium, ammonium, and sodium hydroxides. Assuming two exchange sites, C~ 
and E~, in each clay, the observed titration curves indicate that hydrogen is more firmly 
bound to the E> sites than to the C> sites in the order bentonite > beidellite > “illite”. 
The E> site, or substrate, in the titration of an H* substituted clay, is the weaker acid, 
that is, it is occupied only after all C~ sites are occupied by other cations, Differences in 
bonding strength of the C” sites for the several clays are very small. 

The exchange constants indicate that potassium and ammonium react in a similar 
manner and replace hydrogen more easily than sodium replaces hydrogen, In addition 
potassium, ammonium, or sodium will replace one another more readily than they will re- 
pla hydroge n in the clays. 


? 


INTRODUCTION 

Marshall and Bergman (1942a, p. 54. 57), in their experimentation with 
the electrochemical properties of mineral membranes, have used the mem- 
branes in suspensions of the hydrogen-saturated form of various clays; the 
Putnam beidellite. a bentonite from Wyoming, an “illite” from Illinois, and a 
kaolinite from South Carolina. They titrated the suspensions with potassium 
hydroxide, measuring pH and K* activity electrometrically at intervals during 
the neutralization. Titration curves were then prepared from the data. Curves 
were made in a similar manner on suspensions of the same types of clay using 
ammonium hydroxide (Marshall and Bergman, 1942b, p. 328-329) and 
sodium hydroxide (Marshall and Krinbill, 1942, p. 1082, 1083). 

Garrels and Christ (1956) have shown that the titration of the acid form 
of two clays—a beidellite and an “illite’—can be explained on the assumption 
that the clays behave as if they were composed of two distinct clay acids. They 
utilize the data and curves from the potassium hydroxide titration of the 
Putnam beidellite (Marshall and Bergman, 1942a, p. 54) to determine the 
exchange constants for the reaction HC + K*+ + OH =KC+ H,O where 
HC represents the acid form of the clay, and KC the potassium clay complex. 
From the law of mass action, using a, as the thermodynamic concentration of 
each species, HC =a Ky. Multiplying by the activity product of 
[KC] ay- K 
[HC] ax. (1) 


was derived. Kxe-ne is the exchange constant. 


water, the exchange equation 


Garrels and Christ used the following stoichiometric titration relationships 
to obtain Kxe-ne: 


KOH = K+ + KC (2) 
HC, HC + KC + C (3) 
K+ + H?+ C- + OH (4) 


where HCo is the original concentration or equivalent exchange capacity of 
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the clay acid and C~ the concentration of dissociated clay complex, H* and 
OH in equation (4) were considered insignificant (<10°°°). Then values 
for equation (1) were calculated from the data of Marshall and Bergman 
(1942a) for ax+, ay-, and KOH. A value for HCy (equation 3) was assigned 
which would give the most nearly constant Kxe-ye. 

A titration curve could then be calculated using the appropriate HCo, 
Kxe-ne and the observed data of Marshall and Bergman (1942a) for ax+ and 
ay+ in the above equations. This calculated curve provided a fairly close fit 
with part of Marshall and Bergman’s observed titration curve but deviated 
considerably from the high pH end. 


9 
6 
02 03 04 05 06 O7 08 o9 
NHgOH (mols per liter) 
Fig. 1. Titration curve for hydrogen beidellite, 8 percent, by weight, suspension. 


ig. 
(Marshall and Bergman, 1942b, p. 328). 


Garrels and Christ (1956) then made the assumption that this calculated 
curve represented neutralization of the first of two distinct clay acids in the 
suspension. Using the same method and equations as before, an exchange 
constant was calculated for the second acid. Equations were then set up as 
follows using KE to represent the second potassium clay complex, E~ the con- 
centration of the second dissociated clay complex, HE the second clay acid, 
and HE» the concentration or equivalent exchange capacity of the second acid. 


] an 
= 10 
[KE] a, Ke 10 
[HE] ax 
KOH KC KE + K 
HC, HC KC+ C 0.0835 
HE, HE KE +E 0.03815 
H+ +K C- + E- + OH- 


Values were calculated for KOH in the high pH range and compared to 
the observed curve of Marshall and Bergman (1942a). The combined calcu- 
lated curves closely followed their observed curve except for a small deviation 
in the lowest pH region. 


- 
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The derived exchange constants were taken to indicate the presence of two 
substrates or exchange positions in the clay mineral structure, which were 
designated tentatively by Garrels and Christ as the C~ and E~ sites. 

Garrels and Christ also made similar calculations, using the data of 
Marshall and Bergman (1942a), for the titration of a 10 percent, by weight, 
suspension of “illite” by KOH. The HCp was 0.02525 with a Kxe-nc of 10° 
and HE, of 0.009 with a of 

This paper is a follow-up of the work of Garrels and Christ (1956) and 
the results were obtained using their method of calculation as described above 
from their paper. 

Coleman and Harward (1953, 1954) have set forth an alternate thesis to 
explain the neutralization curves for some clays. They propose that clays 
treated with hydrogen exchange resins or leached rapidly with LV HCl are 
hydrogen clays, whereas electrodialized clays or those treated with dilute acid 
are hydrogen-aluminum clays. Therefore when titrating a hydrogen-aluminum 
clay the hydrogen is neutralized first, followed by a buffering action at pH 6 
to 7 at which point the titration curve levels off. A second inflection than occurs 
in the curve when the exchange of aluminum takes place above pH 7. This in- 
fers that two separate exchange reactions take place with two different cations 
rather than two reactions at different energy levels with one cation, 

The clays used as examples for the calculations in the present paper were 
prepared by electrodialysis, therefore the possibility of aluminum entering the 
reaction cannot be ignored. However most of the titration curves shown here 
tend to flatten out below a pH of 6.4 unlike those of Coleman and Harward 
(1953. p- 6046). 

EXCHANGE CONSTANTS 

A titration curve of an 8 percent, by weight, suspension of Putnam 
beidellite with NH,OH was prepared, using the data of Marshall and Bergman 
(1942b, p. 328), as shown in figure 1. Calculated values of NH,OH, assuming 
one exchange reaction, are shown in figure 2. A large deviation in the high 
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NH4OH (mols per liter) 
Fig. 2. Calculated curve (triangles) versus observed curve (circles). Titration for 
8 percent, by weight, suspension hydrogen beidellite assuming one exchange reaction: 
NH.C] an 
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pH range is similar to that obtained by Garrels and Christ (1956, p. 374) 
for KOH, when they assumed one exchange reaction. When two exchange re- 
actions were assumed the calculated titration curve for NH,OH corresponded 
closely to the observed curve (fig. 3). Deviation at the lowest pH values is 
probably related to the effect of hydrogen saturated clay complexes on the 
glass electrode, as suggested by Garrels and Christ (1956, p. 375). However, 
Ganguly and Mukherjee (1951, p. 1429) have suggested that micas, “‘illites”, 
and montmorillonites may have three exchange sites: outer plane surfaces, 
edges, and interlayer. Perhaps then the deviation at very low pH (about 5 or 
lower) may be due to the early release of hydrogen from very loosely held 
outer plane surfaces at a third substrate. No, attempt was made by me to cal- 
culate a third exchange constant for this site. 


100 
s 
90 
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Fig. 3. Calculated curve (triangles) versus observed curve (circles). Titration for 
8 percent by weight, suspension of hydrogen beidellite assuming two exchange reactions: 


(NH,C] an’ [NH,E] an* 


[HC] axa’ 10 and THE) 


Another alternative to explain the deviation may be the affect of low pH 
on aluminum. From pH 4.0 to pH 8.5, aluminum has extremely low solubility. 
Below pH 4.0, aluminum is relatively soluble (Magistad, 1925). Possibly some 
aluminum exchange takes place at this time in addition to the normal hydro- 
gen exchange. 

Exchange constants were determined and titration curves were calculated 
in the same manner as above for other cation-clay combinations as follows: 

1) Sodium—7 percent. by weight, suspension of Putnam beidellite (fig. 

4) from data by Marshall and Krinbill (1942, p. 1081). 
2) Ammonium—10 percent, by weight, suspension of “illite” (fig. 5) 
from data by Marshall and Bergman, (1942b, p. 328). 
3) Sodium—5 percent, by weight, suspension of “illite” (fig. 6) from 
data by Marshall and Krinbill (1942, p. 1083). 
Potassium—3.277 percent, by weight, suspension of a Wyoming 
bentonite (fig. 7) from data by Marshall and Bergman, (1942a, 


p. 55). 
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5) Ammonium—] percent, by weight, suspension of a Wyoming bentonite 
(fig. 8) from data by Marshall and Bergman (1942b, p. 328). 

6) Sodium—2.8 percent, by weight, suspension of a Wyoming bentonite 
(fig. 9) from data by Marshall and Krinbill (1942, p. 1081). 


ike) 


60 + 


02 03 04 06 08 
NoOH (mols per liter) 
Fig. 4. Calculated curve (triangles) versus observed curve (circles). Titration for 
reent, by weight, suspension of hydrogen beidellite (Marshall and Krinbill, 1942) 
[NaC] an” 2 [NaE] an” 
10 and = 
LHC) [HE] aw,” 


issuming two exchange reactions: i07* 
Curves and exchange constants for potassium-beidellite and for potassium- 
“illite” were the same as those calculated by Garrels and Christ (1956, p. 375, 


AO 
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Fig. 5. Calculated curve (triangles) versus observed curve (circles). Titration for 
10 percent by weight, suspension of “illite” (Marshall and Bergman, 1942b) assuming 


[NH,.C] an* [NH,E] an* 


two exchange reactions: - 10°°* and 
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fig. 3, and p. 376. fig. 4. respet tively). Table 1 lists the exchange constants for 
both C~ and E- sites, the clay acid concentrations for both sites. and the 


ratio of the two concentrations for Putnam beidellite. “illite” and a W yoming 
bentonite as titrated 


against potassium hydroxide, ammonium hydroxide. and 
sodium hydroxide. 
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0.00798 


liter) 


acid concentration 
per 


0.083 
0.0116 
0.0301 


0.0064 


constant 


~ 

=f 


hange 
site] 


(log k) 


low pH [¢ 


Bentonite 


\ 
Ns a NN a 
4 
= 
= 
5 os 
i in N= a 
= 


740 Paul D. Blackmon—wNeutralization Curves and the Formulation 


were determined for the exchange or conversion of HC to HE and the cor- 
responding relation for NH,C to NH,E in the Putnam beidellite. Similar equa- 
tions were set up with the other cations and clay minerals. The constants log 


K, and log K. derived are listed in table 2. 


TABLE 2 
Constants for Exchange or Conversion of HC to HE 
and Corresponding Conversion of Titration Cations 


log K log Kz 
: Conversion or exchange Cation 
Conversion of HC to HE of cations from xC to xE 


Potassium 


Beidellite 


Ammonium 


Sodium 


Potassium 
Ammonium 


Sodium 


Potassiun 
Ammonium 


Sodium 


DISCUSSION 

With the exception of the ammonium cation titration against bentonite. 
the exchange reaction data in tables 1 and 2 show correlating relationships 
among the three cations and the three clay minerals.’ Despite the wide range 
of clay acid concentrations (HCo) for the C~ site, which in effect is the ex- 
change capacity for that site, the constants (log k) are consistent for all three 
clay minerals. They range from —2.4 to —2.6 for the ammonium and potassium 
cations and from —2.8 to —3.2 for the sodium cation. The slightly smaller con- 
stant for sodium with each clay suggests that sodium may not replace the 
hydrogen in the clays as easily as ammonium or potassium do. This seems to 
be true for both C~ and E™~ sites. The constants in table 2 also suggest that 
sodium may be more difficult to replace in the clays than either ammonium ot 
potassium but it is not held as tightly as the hydrogen in either site. 

The clay acid concentration (HE») for the E~ site ranges from about 
1/3 to 1/2 that of the C~ site (table 1). However, there appears to be no con- 
sistent pattern in the ratios of the exchange capacities for these two sites with 
the different cations. Variations in the weight percent suspensions of thes: 
* The ammonium-bentonite titration curves as taken from Marshall and Bergman's data 
(1942b, p. 328, 330) are quite different in their general shape from any of the othe: 
cation-clay curves. The calculated curve closely matches the observed curve (fig. 8) but 


none of the constants derived seem to correlate with any of the other calculated catior 
clay constants. I cannot explain the reason for this anomaly. 
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clay minerals with each cation did not cause much change in the constants or 
the ratios of HCy to HEo. However, total exchange capacities are affected by 
other physical aspects of the clay such as particle size and isomorphous sub- 
stitution (Brindley and MacEwan, 1953). Possibly the ratios of HCo to HEo 
are also controlled by factors such as these. 
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Fig. 8. Calculated curve (triangles) versus observed curve (circles). Titration for 
1 percent, by weight, suspension of bentonite (Marshall and Bergman, 1942b) assuming 
[NHC] an” 10 Ls. an” 10 5.4 
[HC] axne ane THE] 


two exchange reactions: 


Apparently the majority of the exchange positions for each clay mineral 
occur in the C~ rather than the E~ sites. The titration curves suggest that the 
E.~ substrate or site in the titration of an H* substituted clay is the weaker 
acid, hence it is occupied only after all C~ sites are occupied by other cations. 

The exchange constants for the E~ sites clearly indicate that the three 
clay minerals react differently with the three cations. All of the log k constants 
for the E~ sites are much smaller than those for the C~ sites, which suggest a 
much tighter relative bonding of hydrogen at the E~ sites, but the E~ site 
constants occur in three definite ranges. For “illite” they range from -5.7 to 
—5.8, for beidellite from —6.7 to —7.4 and for bentonite from -8.0 to -8.3. Ap- 
parently the hydrogen is more tightly bonded in the E~ sites in bentonite than 
in beidellite with “illite” the most easily replaceable of the three, Some indica- 
tion of this relationship may be noted in the C~ site log k constants for potas- 
sium and ammonium cations where the constants are —2.6 for bentonite, —2.5 
for beidellite, and —2.4 for “illite”, however, the difference is very small and 
may not be significant. As noted before, the ammonium and potassium cations 


seem to have similar reactions whereas the sodium cation appears to replace 
the hydrogen of each clay with more difficulty. This follows the order of sym- 
metry values of Ganguly and Mukherjee (1951, p. 1433-1434). 


The constants for the exchange or conversion of HC to HE (log K,) and 
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for the corresponding conversion of the three cations from xC to xE (log kK. ) 
in the three clay minerals are shown in table 2. The positive constants for k, 
indicate that the hydrogen “prefers” or is more strongly attracted to the E 
than to the C~ substrate. Negative constants for K, indicate a cation “prefer- 
ence” for the C~ sites. The bonding strength difference, however, apparently 
is not as great between the C~ and E> sites for the exchange cations as it is 
for the H+ as shown by the small constants for K,. Garrels and Christ (1956, 
p. 377), in their calculations for potassium-beidellite obtained constants of 
+ 4.0 for K, and zero for K.. Their figures suggest that there is no K* “pref- 
erence’ between the C~ and E~ sites. Both sets of figures indicate, however, 
that the difference in bonding strength between C~ and E~ sites for potassium, 
ammonium and sodium is low. 

It would seem then that total replacement of hydrogen in these three clays 
by potassium, ammonium or sodium would be more difficult than replacement 
of one of the three latter ions by another of them. Further work with divalent 
cations such as calcium and magnesium is necessary to complete the exchange 
reaction picture with these three clays and with other minerals having ex- 
changeable ions in their crystal structures. 


SUMMARY 


Exchange constants have been derived for three ¢ lay minerals. a beidellite. 
an “illite”, and a bentonite. The constants were used to interpret titration 
curves for the three clays using potassium, ammonium, and sodium hydroxides. 


Two exchange sites in each clay. designated C~ and E~ sites are assumed. 
The “exchange capacity” of the C~ sites is from two to three times greater 
than that of the E~ sites. However, hydrogen in the three clays is more firml 
bound to the E~ sites than to the C~ sites in the order bentonite > beidellite 


“illite”. Differences in bonding strength of cations to the Co sites are very 
small among the three clays. 
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Fig. 9. Cal ulated curve (triangles) versus observed curve (circles). Titration for 
2.8 percent, by weight, suspension of bentonite (Marshall and Krinbill, 1942) assuming 
[NaF] 

[HC] a LHE] as,” 


two exchange reactions 
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Of the three cations potassium, ammonium, and sodium, the sodium seems 
to replace the hydrogen with more difficulty than either potassium or am- 
monium, The latter two react similarly. Hydrogen appears to have a “pref- 
erence” or greater bonding strength toward the E~ sites whereas the other 
three cations “prefer” the C~ site slightly more than the E~ site. Probably 
potassium, ammonium or sodium would replace each other more readily in 
these clay minerals than any of them would replace hydrogen. 
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PROPERTIES OF WATER. PART IV. 
PRESSURE-VOLUME-TEMPERATURE RELATIONS OF 
WATER IN THE RANGE 100-400 C AND 100-1400 BARS* 
WILLIAM T. HOLSER and GEORGE C. KENNEDY 


ABSTRACT. 


New precision measurements of the specific volume of water in the range 
100-400 


and 100-1400 bars are generally close to earlier published measurements, Re- 
cent measurements in the critical region by Kirillin, Rumiantsev, and Zubarey are used to 
complete a table of specific volumes for this range. 


INTRODUCTION 

The specific volume of water was measured in the range 100-400 C, 100- 
1400 bars during 1957, as part of the general program of precision measure- 
ments of the pressure-volume-temperature relations of water. Measurements of 
apparently high precision have long been available from the work of Smith 
and Keyes (1934) for the liquid water region to about 350 bars. Kirillin, 
Rumiantsev, and Zubarev (1956) provide a large number of recent measure- 
ments near the critical temperature, extending for some temperatures to 900 
bars. The previous measurements of Kennedy (1950), at temperatures above 
200 C, extended to 2500 bars. Kennedy accepted the data of Smith and Keyes, 
which were in substantial agreements with his own. It was not expected that 
the measurements in this lower temperature region would be subject to ap- 
preciable errors by reaction with the bomb wall as was encountered at high 
temperature and low pressure (Kennedy, 1957). However, plots of isothermal 
compressibility and isobaric thermal expansion indicated some inconsistency 
and lack of continuity in Kennedy’s table; these seem to be most pronounced 
at about 400 bars, where these measurements replace those of Smith and Keyes. 

Precision equipment assembled for the purpose of measurements in a 
higher temperature range (Holser and Kennedy, in press) made it possible to 
check the earlier measurements to 1400 bars, and to extend them down to 
100 C to meet the new measurements in the range 0-100 C (Kennedy, Knight, 
and Holser, 1958). The present measurements are principally concerned with 
the liquid water region. The numerous data by Kirillin were already available 
for the critical region, and an improvement or correction of these values would 
have required much more precise control and measurement of temperature 
than was immediately available. Therefore the Kirillin data have been used 
extensively in completing a smooth table through the critical region. 


APPARATUS 
A schematic diagram of the apparatus used is shown in figure 1. The ar- 
rangement is similar to that used earlier (Kennedy, 1954, p. 230), with the 
addition of a dead-weight piston gauge in the oil system for precise pressure 
measurement. 
Piston gauge corrections for temperature and pressure were made as 
described by Johnson and Newhall (1953, p. 303). The sensitivity of the gauge 
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was about .05 percent. In order to effectively use this sensitivity with the capil- 
lary lines used for connections, it was necessary to use a very low viscosity 
transmission fluid, obtained by mixing equal proportions of kerosene and a 
light oil (Gargoyle Vacuoline type E). The Bourdon gauges are 15-inch Heise 
vauges with a range of 1400 bars. 


FURNACE 
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Fig. 1. Schematic diagram of apparatus. 


The volumometers and mercury U-block have been previously described. 
[he volumometer was calibrated by weighing water ejected at various pres- 
sures, with a result precise to 0.03 percent. The bombs used in these experi- 
ments were similar to those previously described, made of Inconel-X, or 100- 
NT2 alloy. The volumes of about 70 and 170 cc were measured by repeated 
filling with water at about 200 C and 200 bars and ejecting and weighing the 
water; the precision was .03 percent. Volumes at other temperatures and pres- 
sures were calculated from the known thermal expansion and elastic properties 
of the materials. Water was supplied from a flask connected with a reflux con- 
denser in such a manner that steam at one bar was always in contact with the 
water surface, in order to keep gases out of solution. 

Temperature was measured with a Pt, Pt-10 percent Rh thermocouple in- 
serted in a well in the bomb. The thermocouple was calibrated at various times 
against two thermocouples that had previously been calibrated at the National 


Bureau of Standards. The accuracy of temperature measurement was about 
0.2 


METHOD 


Most of the measurements of specific volume were made along isotherms, 
by measuring the compressibility with respect to a low-pressure (300 bar) 
base. The base line was established by measuring the isobaric thermal expan- 
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sion at 300 bars, starting with the previously measured value at 100 C. 300 
bars (Kennedy, Knight, and Holser, 1958), and supplemented by the earliet 
measurements of Smith and Keyes (1934) and Kirillin, Rumiantsev, and 
Zubarev (1956). 

A typical isotherm measurement started with loading of distilled. de- 
areated water into the bomb at the highest pressure, along with a little CuO 
to eliminate errors from hydrogen generated by reaction with the bomb wall. 
Valve 6 was then cut off, so that if the mercury remained at a constant level. 
the mass of water in the bomb was very nearly constant. A correction for th 
change of mass of water in the bomb, due to changes of compression of watet 
in the capillary line, was usually insignificant. The bomb with its contained 
water was brought to thermal equilibrium close to the even “isothermal” tem- 
perature, After a preliminary determination of the pressure using the Bourdon 
gauge, the final measurement was made with valves 1, 2, and 4 open to the 
dead-weight piston eauge. which was balanced for no motion of the mercury in 
the U-tube. Valve 4 was then closed and valve 6 opened, and the volumomete: 
backed off, removing sufficient water to reduce the pressure by 50 to 100 bars. 
The volumometer was read at a calibrated pressure before and after the re- 
moval, with valve 6 closed. After thermal equilibrium had been re-established. 
the new pressure was determined, and the process repeated down to a pressure 
of about 100 bars. Small departures from the even isotherm temperature were 
corrected by approximate values of the thermal expansion of water and of the 
bomb material. 

A slightly different arrangement was used for measurements of the iso- 
baric thermal expansion. The volumometer was used to measure the mass of 
water that had to be removed from the bomb, in order to keep the pressure 
constant as the temperature was raised to successively higher values. A very 
long time would have been required at each temperature, to obtain constant 
temperature over the 15-30 minutes required for a reading of the piston gauge. 
Therefore the pressures were read with the precision Heise Bourdon gauge. 
which in turn was calibrated against the piston gauge before and after the 
runs. Readings were made with the temperature stepping both upward and 
downward. 

An independent check on the measurements of temperature and pressure 
was given by measurements along the saturation line. The saturation pressure 
at a given temperature is known to high precision from the work of Osborne 
and Meyers (1934) ; at 365 C the measured pressure of about 200 bars match- 
ed their value within the 1 part in a 1000 corresponding to the temperature 
uncertainty of 0.1 C. A further check on the system of measurement was pro- 
vided by two isotherms run at 100 C. These measurements correlated closely 


with those obtained previously on similar equipment, but by another operator 
in another laboratory (Kennedy. Knight, and Holser, 1958). 


LOW-PRESSURE BASE LINE 
Originally, we contemplated measurement of only the isothermal com- 
pressibility of water to check Kennedy’s earlier determinations. The specifi 
volumes were to be calculated from the base values at low pressures established 
by Smith and Keyes (1934). It became apparent, however, that the Smith and 
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Keves data should be checked for the following reasons: 
(1) Specific volumes given by formulae 1-6 of Smith and Keyes (1934, 

p. 287-289) depart by as much as 0.003 cc/g from their measured 
values, as shown by their figure 1 and table II; this deviation is much 
greater than the internal consistency of the individual isotherms. The 
formulae were used in making their final tabulation, and the latter 
has been copied into many “steam tables” published since. 
Osborne, Stimson, and Ginnings (1937, p. 427-430) found that above 
330 C, the calorimetric determinations of the saturation volume of 
liquid water differed from those obtained by extrapolation according 
to the Smith and Keyes formulae, although the experimental data 
were shown to be reasonably consistent. 
In the range 0-100 C Kennedy, Knight, and Holser (1958) also noted 
the disagreements of the formulae with previous measurements at one 
bar and with their own measurements at high pressures. 
Recent measurements of Kizillin, Rumiantsev, and Zubarev (1956) 
did not seem concordant with the measurements of Smith and Keyes 
in the range 340-380 C. 

[he question therefore arose as to whether interpolations among the measured 

values of Smith and Keyes are the better values for other temperatures and 

pressures, or whether their formulae provided a real smoothing among the 

isotherms. 

Measurements were made on two runs, at intervals of 5-20 C. The meas- 
urements were continued to 405 C, although the rapid change of thermal ex- 
pansion above 360 C reduced the precision considerably. 

Because of the troubles previously experienced with analytical formula- 
tions, the data were reduced by a graphical smoothing of Av/ AT. The results of 
this smoothing are presented in table 1, along with the derived values for ‘iso- 
haric thermal expansion. Deviations of the experimental points from the smooth 
curve are shown in figure 2; for the most part they are less than deviations 
that might be due to a temperature uncertainty of + 0.1 C. At high tempera- 
tures the deviations of volume caused by a small temperature change, increase 
so fast, that it was found most useful to plot the deviations relative to the 
thermal expansion itself. 

The figure also shows deviations of the values of other investigators from 
our smooth curve. The “measured” values of others represent an interpolation 
(or extrapolation) to 300 bars from nearby actual measurements. To the neces- 
sary or accidental extent that these interpolations may be inaccurate, as they 
certainly are in the critical region, such a comparison does not do justice to 
the data. However, some of the deviations are so large that they cannot be 
ascribed to errors of this sort. 

In general, our smoothed curve agrees more closely with the measurements 
of Smith and Keyes, than it does with their formulation. An exteption occurs 
in the range 300-330 C, where the Smith and Keyes formulation lies closer to 
our curve. One other exception is evident in the Smith and Keyes measurement 
at 200 C, which seems out of line with their other measurements and differs 
appreciably from their smoothed value and from our experimental and ob- 
served values. 
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Taste | 
Specific Volume and Isobaric Thermal Expansion of Water at 300 Bars 
Smoothed Measurements of Holser and Kennedy 


I ov/at 
deg ( ce/g (ce/g-deg C)-10* 
100 1.0290 7.3 
110 1.0366 7.8 
120 1.0447 8.4 
130 1.0533 8.9 
140 1.0626 9.5 
150 1.0724 10.1 
160 1.0828 10.7 
170 1.0938 11.4 
180 1.1055 12.1 
190 1.1179 12.8 
200 1.1310 13.6 
210 1.1449 14.5 
220 1.1597 15.4 
230 1.1756 16.5 
240 1.1927 17.8 
250 1.2111 19.2 
260 1.2310 20.8 
270 1.2526 22.6 
280 1.2762 24.7 
290 1.3023 27. 
300 1.3313 30.5 
305 1.3471 32.5 
310 1.3638 34.5 
315 1.3815 36.9 
320 1.4004 39.4 
325 1.4208 42.3 
330 1.4428 45.5 
335 1.4667 49.6 
340 1.4928 54.5 
345 1.5214 59.6 
35 1.5528 65.7 
355 1.5877 66 
360 1.6273 75. 
365 1.673 99 
370 727 119 
375 1.792 147 
380 1.875 188 
385 1,983 254 
390 2.135 370 
395 2,381 680. 
400 2.815 1010. 


It is not claimed that our measurements are definitive. In the critical re- 
gion particularly the equipment used by the Keyes and Kirillin groups, in- 
volving platinum resistance thermometers and large thermal baths, could by 
expected to afford much more accurate data than our equipment, which was 
essentially designed for operation at higher temperatures, Our new measure- 
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ments may throw some light on the disagreements among previously published 
values. In any event their use as a base for our tables makes our tabulation 
completely independent of all but the well established low temperature, 1 bar 
values used as a base by Kennedy, Knight, and Holser (1958). 

For the isotherms tabulated (table 2) at 380 and 400 C, a lack of any 
nearly data from our measurements (except at 300 bars), has led us to rely 
on Kirillin’s extensive data in deriving these values. Any of the rest of our 
tabulation can be adjusted easily by anyone who prefers different base values 
at 300 bars. 


10 
32 
€ 

o ° 
> 
fe 
+ 
z 4 
-20l— | 

50 100 150 200 250 300 350 


TEMPERATURE °C 


2 
> 
> 
Y 
3 $ 
w ~ 
44 
2 -1.0 
4 4 
i 
3x00 325 350 375 400 425 


TEMPERATURE °C 


HOLSER AND KENNEDY (1958) 
SMITH AND KEYES (1934) 
KEYES, SMITH AND GERRY(I936) 
KIRILLIN, RUMIANTZEV AND ZUBAREV(I956) 
NIEUWENBURG AND BLUMENDAL (i932) 
—— SMITH AND KEYES(I934)(EQUATION AND TABLES) 
veviation TO £0.1°C 


Fig. 2. Differences of specific volume at 300 bars—smoothed values of Holser and 
Kennedy (table 1) less other values. 
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ISOTHERMS 

lsotherms were determined in the above-described manner at 140 (2), 
180, 220, 260, and 298 C. An isotherm at 355 and two at 420 (described in 
Holser and Kennedy, in press) were run all the way down to low pressure where 
the specific volume of the small amount of remaining steam was well known 
from previous measurements (Kennedy. 1957). Therefore these higher tem- 
perature isotherms were more or less independent of the 300 bar base line. 

It was desirable to smooth all isotherms by graphing the data in a man- 
ner that would emphasize its small irregularities. Therefore the method of item 
differences was employed, by which, in effect, the slopes given by the differ- 
ences between succeeding determinations are smoothed, and the smooth values 
integrated back. This procedure is facilitated if the differential (or difference) 
function that is graphed lies on approximately a straight line. If the specific 
volume of water on an isotherm were represented by equations of either the 
Tumlirz (Eckart, 1957) or van der Waals types, 


(p + (w— vo) K (or — RT) (1) 
then the function (-dp/dv)"? should plot as a straight line against pressure: 
dp ‘ 

K p (2) 


On the other hand, if the data for volume is more precisely represented by an 
equation like that of Tait (Hirschfelder, Curtiss, and Bird, 1954, p. 261): 


\ B+ p 
A In (3) 
then the slope (dp dv), rather than its square root, plots as a straight line: 
dp 
Ay + B p (4) 
dy 


Several comparisons in the range 300-360 C indicated that (dp/dv)'? was less 
straight than (dp dv), so the latter function was used in smoothing all the 
isotherms. These curves are slightly concave toward the p axis at high tem- 
perature, but below 300 C no curvature could be detected within the precision 
of the data. A detailed comparison of the various equations of state for liquid 
water will be made in a subsequent publication of this series. 

Most measurements were within + 0.0002 cc/g of the smoothed iso- 
therms, except at 355 C where deviations were as high as + 0.0010 cc/g. The 
individual original measurements are listed and compared with the final tabu- 
lated values (see below) and deposited for reference e.' 

The isotherms in the range 100-360 C were calculated on the basis of the 
300 bar values that are described above and listed in table 1. A change in this 
hase value would affect the calculated values at all other pressures by the same 
amount in cc/g, within the limits of precision of the measurements. 


Available as document no. 5756, American Documentation Institute, Auxiliary Publica- 
tions Project, Photoduplication Service, Library of Congress, Washington 25, D, C. A copy 
may be secured by citing the Document number and by remitting $1.25 for photoprints, 
or $1.25 for 35 mm. microfilm. Advace payment is required. Make checks or money orders 
payable to: Chief, Photoduplication Service, Library of Congress. 
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TABULATION OF RESULTS 

In order to facilitate interpolation between isotherms, and to bring out 
any errors or any necessity for smoothing, the approximate thermal expansions, 
Av/AT, were plotted against temperature for each isobar. In only one place, 
near the upper end of the 180 C isotherm, was any smoothing between iso- 
therms indicated by these plots. Interpolations could be made with certainty 
to derive a table extending to 360 C. With the help of the data of Kirillin, 
Rumiantsev and Zubarev, the best values for the isotherms at 380 and 400 C 
could be deduced. The final result is shown in table 2. 


COMPARISON WITH THE RESULTS OF OTHER INVESTIGATORS 

The results along the 300 bar base line are compared with other measured 
values, and shown in figure 2. In this connection the values we list below 300 
bars cannot be given much weight, inasmuch as we did not measure enough 
points at low pressure. In fact, it was necessary to depend upon the data of 
Smith and Keyes to indicate the detailed trend at low pressure above 300 C. At 
lower temperatures the values listed represent essentially a linear extrapola- 
tion of the Ap/Av plot. At all temperatures the values of Ap/Av derived from 
the measurements of Smith and Keyes are slightly but consistently above ours, 
and those of Kirillin, Rumiantsev, and Zubarev are too irregular to make any 
valid comparison. In the range 100-300 C our measurements of the volume at 
300 bars (table 1) generally lie a little above those tabulated by Smith and 
Keyes, and therefore values of the saturation specific volume derived by our 
extrapolation would be even higher than those of Smith and Keyes. However. 
our data are still consistent with those derived by calorimetry. The calorimetric 
measurements by Osborne, Stimson, and Fiock (1930) and of Osborne, 
Stimson, and Ginnings (1937) give a function B, from which it is possible to 
calculate the specific volume of liquid water: 


B 


T(dp aT) xc (5) 


Vs 


using only the well-established slope of the saturation line. These authors ac- 
cepted the formulated and tabulated values of Smith and Keyes in making 
their final correlation (Osborne, Stimson, and Ginnings, 1939), but an evalua- 
tion of their original data gives specific volumes generally above those tabu- 
lated by Smith and Keyes. The difference amounts to about 1 or 2 standard 
deviations of the calorimetric data. Although this cannot be considered parti- 
cularly significant, it is in the right direction to conform with our extrapola- 
tion. 

At the lower end of the temperature range studied, the data are consistent 
in magnitude and trend with recently published data for the range 0-100 C 
(Kennedy, Knight, Holser, 1958). 

The values reported earlier by Kennedy (1950) for temperatures up to 
300 C are closely confirmed by the present measurements, even at 1400 bars. 
At 400 C the new values are somewhat lower at low pressures and higher at 
high pressures, in line with a trend measured on the isotherms at 420 C and 
above (Holser and Kennedy, 1958). 
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Although the specific volumes measured by Kiyama, Kinoshita and 
Kitahara (1955) at 300 bars (380 and 400 C) are so much greater than our 
smoothed values that they cannot be plotted on figure 2, agreement is better 
with us at high pressure. The single isotherm by Amagat (1893) at 198 C 
differs greatly from all other data, although his measurements below 100 C 
were generally concordant with ours. The measurements of Tamman and 
Ruhenbeck (1932) are very much smaller than ours at 200 C, and very much 
larger at high temperatures. It would appear unwise to give their data any 
weight. 

ACKNOWLEDGMENTS 


This work was supported by the Office of Naval Research under contract 
Nonr 233(27). Parts of the dead-weight piston gauge were kindly loaned by 
F. G. Keyes, and much of the rest of the equipment was made by Lester 
Emerine. Donna Price and G. D. Robinson have kindly read and criticized 
the manuscript; their aid is gratefully acknowledged. 


REFERENCES 


\magat, E. N., 1893, Elasticite et dilatibilite des fluids: Ann, de chim, et phys., ser. 6, 
v. 29, p. 68-136; 505-574. 


Eckart, Carl, 1958, The properties of water. Part Il. The equation of state of water and 
sea water at low temperatures and pressures: Am. Jour. Sct., v. 256, p. 225-240. 


Hirschfelder, J. O., Curtiss, C. F., and Bird, R. B., 1954, Molecular theory of gases and 
liquids: New York, John Wiley and Sons. 


Holser, W. T., and Kennedy, G. C., in press, Properties of water. Part V. Pressure-volume- 
temperature relations of water in the range 400-1000 C and 100-1400 bars: Am. Jour. 


Johnson, D. P., and Newhall, D. H., 1953, The piston gauge as a precise pressure-measur- 
ing instrument: Am. Soc. Mech. Engrs. Trans., v. 75, p. 301-310. 


Kennedy, G, C., 1950, Pressure-volume-temperature relations in water at elevated tempera- 
tures and pressures: Am. Jour. Sct., v. 248, p. 540-564. 


1954, Pressure-volume-temperature relations in CO: at elevated temperatures 
and pressures: Am. Jour. Sct, v. 252, p. 225-241 


, 1957, Properties of water. Part I, Pressure-volume-temperature relations in 
steam to 1000 C and 100 bars pressure: Am. Jour. Sct., v. 255, p. 724-730. 


Kennedy, G, C., Knight, W. L., and Holser, W. T., Properties of water. Part IIT. Specific 
volume of liquid water to 100 C and 1400 bars: Am. Jour. Sct., v. 256, p. 590-595. 


Keyes, F. G., 1933, Methods and procedures used in the Massachusetts Institute of Tech- 
nology program of investigation of the pressures and volumes of water to 460 C: 
Amer. Acad. Arts and Sci. Proc., v. 68, p. 505-564. 


Keyes, F. G., Smith, L. B., and Gerry, H T., 1936, The specific volume of steam in the 
saturated and superheated condition: Amer. Acad, Arts and Sci. Proc., v, 70, p. 319- 
364. 


Kirillin, V. A.. Rumiantsev, L. 1, and Zubarev, U. N., 1956, Experimental investigation of 
the specific volumes of water and steam at high parameters: Fifth International Steam 
Conference, London, Collected Reports of the U.S.S.R. 

Kiyama, Ryo, Kinoshita, Hideo, and Kitahara, Shigeto, 1955, State diagram of steam: Rev. 
Phys. Chem. Japan, v, 25, p, 21-24. 

Meyers, C. H., and Jessup, R. H., 1931, A multiple manometer and piston gauge for pre- 
cision measurements: U. S. Nat, Bur. Stds. Jour. Res., v. 6, p. 1061-1090. 


Nieuwenburg, C. J. von, and Blumendal, H. B., 1932, The isotherms of water from 35- 
480 C and for pressures up to 600 kg./cem.: Rec. des travaux chim, des Pays-Bas, v. 
51, p. 707-714. 


754 illiam Holse r and George Kennedy 


Osborne, N. S., and Meyers, C. H., 1934, A formula and tables for the pressure of satu 


rated water vapor in the range 0 to 374 C: U. S. Nat. Bur. Stds. Jour. Res., v. 13, 


p. 1-20 
Osborne, N. S., Stimson, H. F., and Fiock, E. F., 1930, A calorimetric determination of 


thermal properties of saturated water and steam from 0 to 270 C: U. S. Nat. Bur 


Stds, Jour, Res., v. 5, p. 411-505 
Osborne, N. S., Stimson, H. F., and Ginnings, D. C., 1937, Calorimetric determination of 
the thermodynamic properties of saturated water in both the liquid and gaseous state 
from 100 to 374 C: U.S. Bur. Stds. Jour. Res., v. 18, p. 389-448. 
1939, Thermal properties of saturated water and steam: U. S. Nat, Bur. Stds 
Jour. Res,, v. 23 p 261-270 
Smith, L. B., and Keyes, F. G., 1934, The volumes of unit mass of liquid water and their 
ssure and temperature: Amer, Acad. Arts and Sci 


correlation as a function of pre 
Proc.. v. 69, p 286-312 


lamman, G., and Ruhenbeck, A., 1932, Die spezifischen Volumen des Wassers zwischet 


20 und 650 ¢ Ann. der Phys., v. 405, p. 63-79. 
W. T. Horser 
CALIFORNIA Researcu Cori 
Box 146 
La Hapra, CALIFORNIA 
G. C. Kennepy 
INstTITUTE OF GEOPHYSICS 
UNIVERSITY OF CALIFORNIA 
Los ANGELES. CALIFORNIA 


Vou. 256 No. 10, 1958 


American Journal of Science 
DECEMBER 1958 


CONTENTS 


Stratigraphy and Radiocarbon Dates at Searles Lake, California 
Richard Foster Flint and W. A. Gale 


Glacial Flutings in Central and Northern Alberta 
C. P. Gravenor and W. A. Meneley 


A Further Attempt to Determine the Vapor Pressure of Brucite 
W.S. Fyfe 


Neutralization Curves and the Formulation of Monovalent Cation Ex- 
change Properties of Clay Minerals ............0.-s0+ -Paul D, Blackmon 


Properties of Water. Part IV. Pressure-Volume-Temperature Relations of 
Water in the Range 100-400 C and 100-1400 Bars 
William T. Holser and George C. Kennedy 


Inpex To Votume 256 


Contents, Votume 256 


715 


729 


733 


744 


